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1. Introduction. 



C " 3 Abstract. A new derivation of the classic asymptotic expansion of the n-th 

^Nj , prime is presented. A fast algorithm for the computation of its terms is also 

given, which will be an improvement of that by Salvy (1994). 

O , ' Realistic bounds for the error with li~^ (n), after having retained the first m 

J^ ^ terms, for 1 < m < 11, are given. Finally, assuming the Riemann Hypothesis, 

^^ , we give estimations of the best possible ra such that, for n > rs,, we have 

Xf^ ' Pn > S3(n) where s^in) is the sum of the first four terms of the asymptotic 

expansion. 

H 

C^ I l-l- Historical note. Chebyshev failed to fully prove the Prime Number Theorem 

(PNT), but he obtained some notable approximations. For example, he proved that 
for every natural number n: if the limit 

<y' lim i^i^(4x)-li(x)) 

~^ , exists, then this limit must be equal to 0. 

.^ ■ The question was decided by de la Vallee Poussin (1899) when he gave his bound 

\f^ , on the error in the PNT: The above limits exist and equal 0. 

In 1894, Pervushin, a priest in Perm, published several formulae obtained em- 
("^ , pirically about prime numberqj. One of these formulae gives the following approx- 

^SJ ' imation to the n-th prime 



— = log n + log log n - 1 + — — 5— . 

n '^ '^ ^ 121ogn 241og^n 

Cesaro then published a note [1, 1894] where he asserts that the true formula is 



H 

•^ ■ P« 1 , 1 1 1 , loglogn-2 

■ - - — = log n + log log n — 1 



n log n 

(loglogn)2-61oglogn + ll , ,, _2 x 

— ^ + o(log n). 

2 log n 

Despite no mention by Cesaro in [1], the editors of his collected works added a 
note to \Xi pointing out that certain formulae quoted by Cesaro, since they followed 
from the results of Chebyshev, were only established under the assumption of the 
existence of the implied limits. It therefore remains unsurprising that Hilbert, in 
the Jahrbucto stated that Cesaro did not prove his formula. 



Ivan Mikheevich Pervushin (1827-1900) (HBan Mhxccbhh IIcpByniHH). No small achieve- 
ment if we note that he had only a table of primes up to 3 000 000. 

Jahrbuch iiber die Fortschritte der Mathematik (1868-1942), a forerunner for the Zentralblatt 
fiir Mathematik, at present digitalized at http://www.emis.de/MATH/JFM/JFM.html 
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2 JUAN ARIAS DE REYNA AND JEREMY TOULISSE 

Landau J7, 1907] several years later was better informed: a formula, like that of 
Cesaro, would imply the PNT, which had yet to be proved at Cesaro's time. How- 
ever, using the results of Chebyshev, Cesaro may claim that if there is some formula 
for Pn correct to the order n(logn)~^, then it must coincide with his formula. 

Cipolla 1902] obtained an infinite asymptotic expansion for p„ and gave a 

recursive formula to compute its terms. He published after the results of de la 

Vallee Poussin but it seems that he was unaware of these results, so that gave 

his proof under the same hypotheses as Cesaro. So uninformed was he that he 

attempted to prove some false formulae of Pervushin already corrected by Torelli 

[22] 

loglogn-1 ,^/loglogn\2 

Pn+i - Pn = log n + log log n H h e> — 

logn V logn / 

with an impeccable proof that if such a formula exists, then it must be this formula. 
(Such a formula would refute the twin prime conjecture, and today the above 
formula is known to be false.) 

In his Handbuch Landau [SI § 57] obtained by means of the procedure of Cesaro, 
some approximative formulae for p„, and explained that the method could give 
further terms. He also mentioned some recursive formulae without giving any clue 
for their derivation. 

We may say that Pervushin was the first to deal with a formula for p„, albeit 
that he gave only the first few terms. Cesaro then proved that in the case such 
a formula exists, it must be one from which he would be able to derive several 
terms. Cipolla found a method to write all the terms of the expansion if there is 
one. Landau saw that the results of de la Vallee Poussin imply that the expansion 
certainly exists. 

The algorithm given by Cipolla is not very convenient for the computation of the 
terms of the expansion. He iteratively computes the derivative of some polynomials 
appearing in the expansion but computes the constant terms as determinants of 
increasing order. Robin [ITl 1988] considers the problem of computing these and 
other similar expansions, leaving the problem of computing the constant terms of 
the polynomials as an open problem. Later Salvy [20l 1994] gives a satisfactory 
algorithm. This algorithm needs O{n7'^y/\ogn) coefficient operations to compute 
all the polynomials up to the n-th polynomial. 

The asymptotic expansion of p„ also plays a role in the study of g{n), which is 
the maximum order of any element in the symmetric group 5„. In fact, log g(n) 



has the same asymptotic expansion as y li (x) |llj . 

There are many results giving true bounds on p„, for example we mention p„ > 
nlogn [18J 1939], and p„ > n(logn + log logn - 1) [41 1999] both for n > 2 (with 
partial results given in [19], [16], [12], [4]). In [5] it is also proved that 

/ log log n — 2 \ 

Pn < n\ log n + log log n — 1 H , n> 688 383. 

V log n / 



1.2. Organization of the paper. In this paper we present a new derivation of 
the asymptotic expansion for p„ and obtain explicit bounds for the error. 

First, it must be said that the asymptotic expansion has, in a certain sense, 
nothing to do with prime numbers: it is an asymptotic expansion of ali(x) := li~ x 
which is the inverse of the usual logarithmic integral function. 
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In Section [3] a proof of the existence of the expansion is given, fohowing the path 
of Cesaro, since it cannot be found elsewhere, ahhough it is frequently claimed it 
can be done. This Section is not needed in the rest of the paper. 

In Section|4l a new formal derivation of the expansion is given. We obtain a new 
algorithm to compute the polynomials (Theorem 14.91) . This is simpler than that 
given by Salvy [20' . Our algorithm allows all the polynomials up to the n-th one to 
be computed in O(n^) coefficient operations (Theorem 14. 11 1) . It must be said that 
these polynomials have 0{'n?) coefficients. 

In Section [SJ independently of Section [31 we prove that the formal expansion 
given in Section |4] is in fact the asymptotic expansion of ali(a;) and gives realistic 
bounds on the error (Theorem 15. 161 and 15. 17|) . 

In Section [6l the results are applied to p„ the n-th prime. Using the results of 
de la Vallee Poussin it can be shown that the asymptotic expansion of ali(ri) is also 
an asymptotic expansion for p„ . 

By assuming the Riemann Hypothesis, we found (Theorem I6.2p that 

\pn- &\\{n)\< -y/n{\ogn)i , n > 11. 

TT 

This bound of p„ is better than all the bounds cited above. 

We end the paper by motivating why the above bounds have not been extended 
to further terms of the asymptotic expansion (Theorem 16. 4|) . 

Notations: With a certain hesitation we have introduced the notation ali(a;) to 
denote the inverse function of li(a;). 

In Section [51 where explicit bounds are sought, it has been useful to denote by 
9 a real or complex number of absolute value \9\ < 1, which will not always be the 
same, and depends on all parameters or variables in the corresponding equation. 

Acknowledgement: The authors would like to thank Jan van de Lune ( 
Hallum, The Netherlands ) for his linguistic assistance in preparing the paper, and 
his interest in our results. 

2. The inverse function of the logarithmic integral. 
Usually li(a;) is defined for real x as the principal value of the integral 

h(x) = p. V. / . 

^ ' Jo log^ 

It may be extended to an analytic function over the region Q ~ C\ (— oo, 1], which 
is the complex plane with a cut along the real axis x < 1. The main branch of 
the logarithm is defined in fl and does not vanish there. Therefore, li(z) may be 
defined in fl by 

(1) ii(^)^ii(2)+ r^, zen 

where we integrate, for example, along the segment from 2 to z. 

For real a; > 1, the function li(a;) is increasing and maps the interval (1,+cx)) 
onto (— oo, +oo), so that we may define the inverse function ali: R — )■ (1, +cx)) by 

(2) h(ah(x)) = X. 
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The function li(2:) is analytic on f2, so that ali{x) is real analytic. It is clear that 
we have the following rules of differentiation 

(3) — \i{x) = , — ah(a;) = logah(a;). 

ax log X ax 

It is well known that the function \i{x) has an asymptotic expansion: 
Theorem 2.1. For each integer N > 

logxV ^log'^a; Mog^+ia;yy' 

This may be proved by repeated integration by parts (see [13l p. 190-192]). 

3. Asymptotic expansion of ali(a;). 
In this section, we prove the following 
Theorem 3.1. For each integer N > 

n—l 

where the P„_i(z) are polynomials of degree < n. 

In the case of A^ = the sum must be understood as equal to 0. 
The theorem says that, for each N, there exists an xn > 1 and a constant Cn 
such that 

ah(e^) ^ .^P„_i(logx) '-^+1 



V^ -r„-Uiogx-j log X 

1-Z^ Z^ ^^^ ^N+l ' {x>xn). 



xe-^ 

n— 1 

In the course of the proof we will make repeated use of the following 
Lemma 3.2. Let f{x) be a function defined on a neighbourhood of x = such that 
(6) f{x)^aix + --- + aNX^ + 0{x^+'^), {x -^ 0) 

where the ak are given constants. Assume that g{x) satisfies 

N n \ 1 N+1 

f7\ ( \ \^ Pni^Ogx) .^f^Og_X\ t ^, \ 

n—l 

where the Pn{z) are polynomials of degree < n. Then there exist polynomials qk{z) 
of degree < k such that 

^ n \ 1 N+1 

(8) f^g(x))^^^—^ + 0^ j, (X^+OO). 

Proof. It is clear that, for each 1 < n < A we have p„(loga:)a;~" = O((loga:/x)"). 
Therefore, g{x) = O{\ogx/x) (this is true even when A = and there is no p„). It 
follows that linij; _>._(_ 00 .9(2;) = and by substitution in ([5]), we obtain 

N -, Af+l 

(9) /(g(x)) = 5]a„g(x)" + 0^°^ "^ 
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By expanding the powers g{x)^ by ([7]) it is easy to obtain an expression of the form 

^ n \ 1 N+1 

(10) gix) ^^ -^ + e)(^___j, (x^+oo) 

fc=l 

where each Pn.k{z) is a polynomial of degree < k. By substituting these values in 
equation ^ and collecting terms with the same power of x, ([5]) is obtained. D 

We will prove Theorem 13.11 by induction. The following theorem yields the first 
step of this induction. 

Theorem 3.3. 

(11) ^^,^aC-^^), {X-.+00). 

a; log a; V logx / 

Proof. From Q with iV = we have "'^\'°g^ = 1 + O(log"^ y) for y ^ oo. Since 
limj;^_i_oo ali(a;) — +oo we may substitute y ~ ali(2:) and obtain 

, , X log ali(a;) 

(12) ° ^ =1 + 



ali(a;) Vlogali(.T) 

By taking logarithms 

log X - log ah(a;) + log log ali(x) = 0( -— — 

Vlogali(x) 

we obtain 

log a; _, , ,^/loglogah(x; 



logali(2;) V logah(a;) 



(13) r^ = i + K- 

logah(a;) V 

and it follows that 

(14) lim ^^ = 1. 

x-)-+oo log ali(x) 

By taking log in P^ 

'log log ah (a;) 



log log X — log log ali(a;) — 0(- 



log ali(x) 



we obtain 


loglogx I \ oi ^ 




loglogali(x) ' "'Vlogali(x) 


so that 




(15) 


log log a; 
Jl™ ^_,_^,^, , =1. 



In view of p4)) and p5)) . we may write (|12l) and (1131) in the form 

^^g^ xlogah(x) _^ ^ gf M, log^ _i I ^>;>^loglog^ 



ali(a;) Vloga;/ logah(a;) V log a; 

and by multiplying these two, we obtain 

xlogx _ _^ ^^loglogx 



ali(a:) V logo; 

from which (jlip can easily be deduced. D 
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Proof of Theorem \3.1\ We proceed by induction. For A^ = 0, our theorem is simply 
Theorem 13.31 with e^ instead of x. 

Hence we assume ([S]) and try to prove the case A'^ + 1. 

Our objective will be obtained by starting from the expansion of li(y). By ([4]) 

By substituting y = ali(e^) and applying (ITi)) we obtain 

(17) e-logali(e-)^ ^^ k\ ^o(^—\ 

^ ' ah(e^) ^ (logah(e^))fc Va;^+2y'- 

From our induction hypothesis, the expansion of logali(x) and (logali(a;))~'^ is 
now sought. 

By taking the log of ^ we obtain 

,og.l.(.') = . + log. + log{l + £ ^l^^fifi + 0(^1;;;^;)}. 

n— 1 



Lemma 13.21 may be applied with 



log(l + X) = X-^ + ^~--- (-1)^+% + 0(X~+i) 



-^ Q„+i(log.x) , ^^/log^+^x^ 



to obtain 

logali(e )=a; + logx + ^ — + 0[^ ^^^^ j. 

n—l 

The reason why we have written Qn+i instead of Qn{x) is revealed below. The 
above may be written as 

logah(e)^x|l + ^ + )^ ,;»+i +^(^lvT^J| 

Tl=l 



(18) logaH(e^) = a;{l+ ^ 



^+1 ^ n,.„^^ .i^„A^+2 



O"(loga;) , ,^/log™+^x 



O 



X" 
n—l 



rN + 2 



where the Q„(z) are polynomials of degree < n. Observe that knowing the ex- 
pansion of ali(e'^) up to (loga;/a;)''^+^ has enabled us to obtain logali(e'^) up to 
(logx/x)^"'"^; this will be of great importance in what follows. 
From (|18p . for all natural numbers n, 





1 1 


{:. 


N+1 

Tr 

fc=l 


9/c(l0g 


^'+c.(i 


og^+^. 




(logaH(e^))" x" 


^N+2 


By 


applying Lemma |3.2| with 




Af + 1 


/ Tl\ 








(1 + =^)-" 


- 1 = 


r=l 


I :) 


x' + 0{x 


-N-2\ 



)}' 
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we obtain 



1 1 / 

lloffalife^)!" ~ x" V ^ 



^1 T/ .n„„^^ ,i^„A^+2 



Ka(loga:) , ^(\og''+'x 



{logali(e^)}" x" V ^ x*^ V x^+'^ 

where the Vn.k{z) are polynomials of degree < k. 

By substituting these values of {logali(e^)}^" in ([TT]), we obtain 

e"logali(e") ^ , , ^' Ukjlogx) / log^+^x x 

ah(e^) ^^ x'= V x^+2 )' 

Hence again from ([T51) with n = 1 

-'V+l TA /I N 1 JV+2 

— T4,fe(log2;) „/log ^ a;^ 



ali(e^ 



i{,,j;Vi^,,(i^)J 



fc=l 



-'V+l TT n \ 1 -/V+2 

{i+i:^^+c>(i2|^)} 



fe=l 



from which we derive that there exist polynomials Wk{z) of degree < k such that 
xe- _ ^+^Wk{\ogx) , ^^/log^+'' 



^ ^ Wk{^ogx) .^/ iog x \ 

^+Z^ X'= V X^+2 y*- 



ah(e-) 



Another application of Lemma 13.21 yields 

ali(e-)^ ^^ T,(logx) / log^+^x 

xe^ ^ xfe V x^+2 

fc=i 

with polynomials Tk{z) of degree < k. Therefore, we have an asymptotic expansion 
of type ([Sj with A^ + 1 instead of N . The usual argument of uniqueness of the 
asymptotic expansion applies here so that Tk{z) = Pk{z) ior 1 <k < N . D 

4. Formal Asymptotic expansion. 

First we give some motivation. We have seen that the asymptotic expansion of 
ali(e^) is 

an(e^) = xe^F(x, log x), where F(x, y) := 1 + ^ ^n 

n—l 

and differentiation yields 

e^ logah(e^) = (e^ + xe'')V + xe^K + e'^Vy. 
Here logali(e^) — log(xe^y(x,logx)) = y + x + logV^, so that 

y + x + logV = V + xV + xVx + Vy 
which we write as 

(20) V ^l+'^--V-V^--Vy + -\ogV. 

XX XX 

This ends our motivation for considering this equation. 
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Consider now the ring A of the formal power series of the type 

y^ Qnjy) 

where the q-niy) are polynomials with complex coefficients of degree less than or 
equal to n. In particular go(y) is a constant. 

It is clear that A, with the obvious operations, is a ring. The elements with 
go = form a maximal ideal /. An element 1 + m with go = 1 is invertible, with 
inverse 1 — u + u^ — ---. It follows that if a ^ /, then a is also invertible. Hence / is 
the unique maximal ideal and A is a local ring. If a € A is a non- vanishing element, 
then there exists a least natural number n with g„(y) ^ 0. We define deg(a) = n 
in this case, with deg(O) — oo. 

As is usual in local rings, (see [9]) we may define a topology induced by the 
norm ||a|| = 2^'^'^^^°-\ which, with the associated metric, induces a complete metric 
space. Indeed A is isomorphic to C[[X, Y]\, by means of the application that sends 
X !->• x~^ , Y i-> yx~^. 

Given a G A with a — Yll^=o •^ , we define two derivates 



= 1 -^ n=l 



X" 



Finally the set U C A oi elements with qo = I form a multiplicative subgroup of 
A* (the group of invertible elements of A). For 1 + u E U, we define 

oo J^ 

k=l 

which is a series that is easily shown to converge since u'' E I^ . 
We are now ready to prove the following 

Theorem 4.1. The equation (|20p has one and only one solution in the ring A. 

Proof. For V^ G [/, we define T{V) as 

T{V) ■.= l + ^--V-V.,^-Vy + - logK 

XX XX 

It is clear that T{V) € U . We may apply Banach's fixed-point theorem. Indeed, we 
have deg(r(y)-T(VF))<l+deg(F-'l^), so that \\T{V)-T{W)\\ <\\\V-W\\. 
By Banach's theorem there is a unique solution to F == T{V). We may obtain 
this solution as the limit of the sequence T"(l). In fact, since deg(T(F) — T(W^)) < 
1 + deg(y — W), in each iteration we obtain one further term of the expansion. In 
this way, it is easy to prove that the solution is 

w- 1 w-2 

V = 1 + ^ + ^-5- + --- 

X X 

However, we are going to find more direct methods to compute the terms of the 
expansion. D 

Definition 4.2. Let V be the unique solution to equation (|20p . Since it is in A, it 
has the form 

(21) n-,j^)-i+f:^^ 

n— 1 
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where for n > 0, Pn-i(y) is a polynomial of degree < n. 

In the following sections, we prove that V yields the asymptotic expansion of 
ali(e^). For this proof the following property is crucial. 

Theorem 4.3. For N >1, let 

(22) w{x,y)=WN{x,y):=l + f2^^^^- 

k=l 

Then 

Proof. By the definition of the P„ we know that deg(T^ — W) > N + 1. Therefore, 
deg(T/ - T{W)) >N + 2. That is 



V-T{W) 



"0(2/) , vo{y) 



X 



N+2 ' rf.N+3 



where wq, vq are polynomials. We also have 



so that 



That is, 



V 


= w + 


00 

>: 


Pn 


-1(2/) 








uoiy) 

^N+2 + 


Voiy) 

xN+3 


+ ■ 


n 


00 

>; 

,=Af+l 


Pn-1 




X" 



W-T{W) 



w-Tiw)= ^^^^^ ■ "^^^ ■ ""^y^ 



j,Af+l ' x^+2 ' x^+3 

for certain polynomials u, v, . . . D 

In the sequel V will denote the unique solution to (PO)) . The element logF 
belongs to A, so that there are polynomials Qn{y) of degree less than or equal to 
n such that 



(24) fog^-E 



Qniy) 



X" 
n=l 



From equation (|20p . we may obtain logV^ in terms of V and its derivatives. It is 
easy to obtain from this expression the following relation 

(25) Q„(2;)=P„(2/)-(n-l)P„_i(y) + P,;_i(y), (n > 1). 

Theorem 4.4. The polynomials Pn{y) that appear in the unique solution (|2ip to 
equation (|20p may be computed by the following recurrence relations: 



Pq ~ y — 1, and for n > 1 

(26) ^ ""' 

i-„ = ni-n-i - n,-i + 

n 



-. n— 1 

p„ = 7ip„_i - p;,_i + - ^ fc{(fc - i)Pfc_i - Pfc - p^^_ jp„_fe_i. 



fc=i 
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Proof. By differentiating ([M]) with respect to a;, we obtain 

/y^ nQn{y) \ (i^Y^ Pn-i{y) \ ^ y^ nPn-i{y) 

n— 1 n— 1 n— 1 

By equating the coefficients of x^"~^, we obtain 

n-l 

(27) nP„_i == nQ„ + ^ kQkPn~k-i, {n > 2). 

fc=i 
Now we substitute the values of the Q„ given in ([25]) 

n-l 

nP„_i = 7iP„ - n{n - l)P„_i + nP,j„i + ^ kQkPn-k-i 

k=l 

so that 

n-l 

nP„ - n^Pn-i - nP;^_, - ^ k{Pk - (fc - l)Pfe-i + P^_i}P„_fc_i. 

fc=i 

D 

From this expression it is very easy to compute the first terms of the expansions 

y-1 y-2 y2_6j/ + ll 2y^ - 21y'^ + 84y - 131 



V = 1 



2x3 6a;4 

3y4 - 462/3 + 294J/2 _ QMy + 1333 



12x5 
y - 1 y^ - 4j/ + 5 2y3 - I5y2 _|_ 42^^ _ 47 



2x2 6x3 

3y4 - 34y3 + i56j/2 _ 366y + 379 



12x4 
Theorem 4.5. M^e /lawe 

(a) For n> I, the degree of Pn is less than or equal to n 

(b) n\Pn{y) has integer coefficients. 

Proof. The equation (j20p may be written 

V-l-^^- (log V-V-xV^~Vy). 

XX 

Since xVic S A, it is clear that 

" Pn{y) 



xV-x-y=-l^Y.-^eA. 



X" 
n—1 



This implies that the degree of P„ is less than or equal to n. 

We prove (b) by induction. The first few P„ satisfy this property. We define 
Pk :— k\Pk so that the recurrence relation (125]) may be written as 



Pn = n^Pn-i - npj^_^ + 



'n — 2 



Hence, by induction, all Pn have integer coefficients. D 
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The most significant contribution by Cipolla is his proof of a recurrence for the 
coefRcients an,k of P„ (see ([30])), which is better than the recurrence given in (f26| . 
We intend to give a shghtly different proof. The result of CipoUa is equivalent to 
the following surprising fact: The solution V of equation (pO)) formally satisfies the 
following linear partial differential equation: 

(28) V^{x-l)Vy-xV^. 

This equation can easily be deduced from the following Theorem. 

Theorem 4.6. For n> 1, we have 

in - i)Pn-i{y) = K-iiy) - Kiy), in > i) 
in~i)Qn-iiy)^Q'n^M-Q'niy). («>2). 

Proof. We will proceed by induction. For n < 3 it can be verified that these 
equalities are satisfied. 

We now assume that (|29p is satisfied for n < N, and we will show that these 
equations are true for n = N + 1. 

By differentiating ([24| with respect to y we get 



/y- Qniy) \ f 1 , y^ Pn-iiy) \ ^ y^ K-iJy) 

Ti— 1 n— 1 n— 1 

so that by equating the coefficients of x^^^^ and of x^^ we obtain 

Af-l N-2 

Qn+1 — Pn ^ /_^ PkQN-k^ Qn — Pn-1 ^ /_^ PkQN~k-l- 

k=0 k=0 

Subtracting these equations we get 

N-2 

Qn+1 ^ Qn = Pn ^ Pn-1 ^ /_^ PkiQN~k ~ QN-k-l) ^ Pn-1 

k=0 

and by the induction hypothesis this is equal to 

N-2 



{N ~ 1)Pn-i +J2Pk-iN-k- l)QN-k-i - Pa 

= -NPn-1 + J2 kQkPN-k^i- 



k=0 

N-1 



fe=l 

By (|27|) this is equal to NP^^i — NQpf so that we obtain 

Qn+1 ~ Qn ~ —NQn. 



This is the second equation of ([29]) for n = iV + 1. In order to achieve the result 
for the first equation, observe that from (P5|) we get 

NQn = NPn - N{N - 1)Pn-i + NPj^_, 
~Q'n = -P'n + iN- 1)P'm-i - P'n-i 
Qn+1 = Pn+1 ^ ^Pn + Pn- 
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By adding these equations we obtain 

- NPn -Pn + P'n+i + ^{^w-i -Pn-{N- 1)^w-i}- 

~ {P'n-i ~P',i-{N- l)P^_,} = NPn -P'm + P'n+i 
which is the first equation of (P^ for n = A^ + 1. D 

We define the coefficients a„ j, imphcitly by 

(30) P„(y) = ^ J, (a„,o;/" - a^ay""' + ' ' ' (-!)"««.„) = 

Analogously, (5„ is of a degree less than or equal to n, and we define the coefficients 
hn,k imphcitly by 

(31) Qn{y)= ^ \ V(-l)H„,fcy"-^ (n>i). 

Remark 4.7. Po{y) has degree 1, which is not given by pop . However, we can 
extend the definition of a{n, k) in such a way that, for n > 1 we have a(n, k) ~ 
for fc < or fc > n. Then a formula such as (|30|) also holds for n == if we add 
up the values from k — —1 to k = n and define a(0, 0) = 1, a(0,— 1) = 1 and 
a(0, fc) = for other values of k 

(30 bis) P„(y)^^^^ ^ (_!)'= a„,y"-'=, (n>0). 

fe=-i 
Note that Qo{y) remains undefined. 

Theorem 4.8. For I < n and < k < n, we have (when defined) 

^(^i ^1) , , ^i(^ ^ 1) , 

(32) an.k = na„^i k-i H ; — an-i,k, bnh — nbn-i,k-i -\ ; — o„_i t. 

n — k n — k 

For 1 < n and < k < n, we have 

(33) On^k = bn,k + {n- k + l)an,k-i- 
For n> 1, we have 

(34) bn.n = nan-l,n-l + /J ( I. ) ^ ^'"''^ "'ri-k-l,n-k-l- 

Proof. ([5^ is obtained by equating the coefficients of y"^-^-^ in the first equation 
in ((29l) . In this way, we obtain 

I'l'^n+k 

(n — 1)! 

^{n- k)— -— — a„_i,fc-i - (n - fc) j a„.fc. 

(n — 1)! n\ 

If n 7^ fc, then the equation for an,k in \Z2\ is obtained. The other equation in bn,k 
is obtained analogously from the second equation in (f29l) . 
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To prove ([33|) . observe that by ([25]), Q„ = Pn — {n ~ l)P„_i + Pn-i, and from 
291) it follows that 



(35) Qn^Pn+P'n, {u > I) . 

Now by equating the coefficient of y"^^^ in both members of this equality we obtain 

Finally (p4)) follows from (l27l) . Recall that — .Sziii and — ^ are respectively the 
values of P„(0), and (5„(0). Hence, by setting y = in ([27)) . we obtain (p4)) through 
multiplication by {n — 1)! and the reordering of the terms. D 



The main problem now is that equations (1321) do not allow us to compute the 
coefficients a„_„. CipoUa gives an algorithm to simultaneously compute the coeffi- 
cients a„.fc and 5„,fc based on Theorem 14.81 In the procedure of Cipolla, these key 
coefficients an^n are recursively computed using all the previous coefficients. We 
prefer a method that computes An := a„_„ and i3„ := ««,«-! separately and then 
compute the remaining coefficients by using p2p . 

Theorem 4.9. In order to compute the numbers a„_fe, we may first compute the 
sequences An '■— On^n o,nd Bn '■— On.n-i by the recursions 

(36) Ao-1, Ai-2, Bo = l, Bi = l, 

(37) Bn = nBn-i+n{n-l)An-i 



n-l ^ „ 

rt — 2 



(38) An = n^An-i + nBn-i- 

n-l 

- (" - 1) E ( r 1 1 ) {^(^ - i)^'«-i - ^'^ + ^^fe-i} A.-fe-i. 

After this one we may obtain a{n,k) := On.k- Setting 

a(0,0) = l, a(0,-l) = l, a(l,0) = l, a(l,l) = 2 
and all other a(0, k) and a(l, k) — 0. Then, for n > 2, put 
a{n,n) = An, 

(39) a{n,k) = na{n-l,k-l) + — -^a(ri-l,fc), (0 < fc < n) 

n — k 

where a(n, fc) = for k < or k > n. 

Finally, we may obtain the b(n, k) :— bn,k from 

(40) 5(71, fc) = a(n, k) - {n - k + l)a{n, k - 1). 

Proof. The constant term of P„ is — ^ and the coefficient of y in P„ is -^, so that 
equation ([57]) follows from the first equation in ([2^ taking it with y = 0. 

In the same way, ([55]) follows from (P5|) . and (15^]) is the first equation in (|32p . 

Equation (HUl) for the b{n,k) follows easily from (|35l) . D 

The array of coefficients a(n, fc) for < n, fc < 7, reads 
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1 























1 


2 




















1 


6 


11 

















2 


21 


84 


131 














6 


92 


588 


1908 


2666 











24 


490 


4380 


22020 


62860 


81534 








20 


3084 


35790 


246480 


1075020 


2823180 


3478014 






720 22428 322224 2838570 16775640 66811920 165838848 196993194 



and the b{n, k) for 1 < n < 7 and < fc < 7 are 



1 


1 




















1 


4 


5 

















2 


15 


42 


47 














6 


68 


312 


732 


758 











24 


370 


2420 


8880 


18820 


18674 








20 


2364 


20370 


103320 


335580 


673140 


654834 






720 17388 187656 1227450 5421360 16485000 32215008 31154346 



Theorem 4.10. (a) The coefficients b(n,k) are integers. 

(b) a{n, fc) > and h{n, k) > 0. 

(c) a{n, fc — 1) < a(n, k) for 1 < k < n. 

(d) Forn> 1, a(n,0) = (n - 1)!. 



Proof, (a) We have proved in Theorem 14.51 that the numbers a{n,k) are integers, 
so that from (HU)) . the coefficients b{n,k) are also integers. 

(b) We proceed by induction on n. Assuming that we have proved that a{m, k) 
and 6(rn, fc) are positive for to < n, it foUows from p2p that a(n, fc) and b{n, k) are 
positive for < fc < n. Then (p4l) imphes that b{n, n) > 0, and (p3| with k — n 
proves that a{n,n) > 0. 

(c) This is a simple consequence of ([55]) . 

(d) The equation follows from (15^ by induction. D 

Theorem 4.11. By means of the rule in Theorem \4.9\ one may compute all coeffi- 
cients Un.k of the polynomials Pn{y) for 1 < n < N in 0{N'^) coefficient operations. 

Proof. We count the operations needed, following the indications in Theorem 14. 9[ 
to compute every a„^fc for < n < iV and < fc < n. 

First we must compute the numbers (™) for < to < iV — 2. Using the scheme 

of the usual triangle, we need to carry out J2k=i ^ additions, which involves {N — 
2)(7V-3)/2 operations. 

The numbers Bn must now be computed for 2 < n < A^ by means of the formula 

Bn^n* {Bn-i + {n - 1) * A„_i). 

Each Bn requires 4 operations, therefore a total of 4(A^— 1) operations are needed. 
We compute the An ioi 2 < n < N using the formula 

An — n * n * An-i + n * Bn-i — [n — 1)* 

* ^ ( , _ j * {fc * (fc - 1) * Ak^i - Ak + k* Bk-i} * An-k-i- 

k=l ^ ^ 
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Hence An requires 7+Y^]^^-^ 8 = 871 — 1 operations. All An together take J2n= 



in- 



1) = AN'^ + 3N — 7 operations. These numbers are the a„^„. The ao^k and ai^k 
require no operations. Finally we compute for < k < n 

an,k = n* {a„_i,fe„i + {n - I) * a„„i^fe/(n - k)}. 
Therefore, a„_fe takes 6 operations. For each n, every Un^k for 1 < k < n takes 
6(n— 1) operations. And each an,k for 2 < n < TV takes X]n=2 6(n— 1) — 3N{N—1). 
The total cost in number of operations is therefore 

^ — ^^ — ^ + 4(iV - 1) + AN'^ + 3N-7 + 3N{N - 1) = 

= ^{15N^+3N-16). 

D 

5. Bounds for the asymptotic expansion. 

5.1. The sequence (a„). First we define a sequence of numbers as the coefficients 
of a formal expansion in A. 

Lemma 5.1. There exists a sequence of integers (a„) such that 

00 , T 00 _, 

(41) log(i-E-) =E--- 

n—l n—1 

The coefficients may be computed by the recursion 

n-l 

(42) ai = 1, On — nl ■ n + y kl On-k- 

fc=i 

Proof. It is clear that u = 1 — J2 nlx^'^ E U C A, so that u^^ E U and logw^^ are 
well defined. To obtain the recursion we differentiate (1411) to obtain 



■^-^ n ■ nl /\-^ a„ \ /_ ^-^ n! 



\l^^n+l)[^ 1^ 



X" 
n—l n—l n—l 



Equation ((42]) is obtained by equating the coefficients of a; " ^. The recurrence 
(|42p proves that a„ is a natural number for each n > 1. D 



The first terms of the sequence (an)5^i are 

1, 5, 25, 137, 841, 5825, 45529, 399713, 3911785, 42302225, . . . 

Lemma 5.2. For each natural number n we have 

(43) a„ < 2n • nl. 

Proof. We may verify this property for ai, 02, 03 and 04 directly. For n > 4 we 
proceed by induction. Assume the inequality for ak with fc < n, so that by (j42p 



rj,' ■ n ^ — ^ 



k /n^ 



nl ■ n ^-^ iji — fc)! • (n — fc) n \k 



< 



k=l 



,fc/ n^ J \n n^ n{n — 1) 

fc=2 ^ ^ ^ ' 
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For n > 4, it is easy to see that this is < 2. 



D 



Lemma 5.3. For each natural number N there is a positive constant cm such that 
(44) a;(i_V:^)>i, x>CN. 



Proof. It is clear that the left-hand side of ((31]) is increasing and tends to +00 when 
X — >■ +00, from which the existence of cjv is clear. 

The value of cjv may be determined as the solution of the equation 

N I 



(45) 



1, 



X > 1. 



In this way we found the following values. 



Cl 


2 


C6 


4.15213 


Cll 


5.61664 


C20 


8.70335 


C2 


2.73205 


cy 


4.43119 


C12 


5.93649 


C30 


12.34925 


C3 


3.20701 


cs 


4.71412 


Cl3 


6.26449 


C40 


16.03475 


C4 


3.56383 


cg 


5.00517 


Ci4 


6.59947 


C50 


19.72833 


C5 


3.86841 


ClO 


5.30597 


Cl5 


6.94035 


C60 


23.42351 



D 



Remark 5.4. Notice that for x > cn the sum in (|44|) is positive and less than 1. 

Proposition 5.5. For each natural number N there exists cIn > such that, for 
a; e C with \x\ > cIn, there exists 6 with \9\ < I such that 

n—1 n—1 



N+lx^+^' 



> dN- 



Proof. By comparing the expansions (j41l) and 

N , , N 



(47) 



it is clear that 



I \ -1 



M-E^)^E 



a„ 1 



E 



&„ 1 



+ (^ + 1)! 



°"^-i^ , so that bj^^i < fljv^i. 



The above expansion is convergent for all sufficiently large \x\, so that 

Ebn ^ _ bN+i 1 , . 

77^ " TV + 1x^+1'^^^^^ 

n=N+l 

where lim^^^oo gN{x) = 1. Hence there exist sufficiently large cIn such that 

\bN+igN{x)\ < ON+l, \x\ > dN- 

This ends the proof of the existence of dN ■ 
We have 



('-(-Ei: 



I \ -1 



n=l 



L. ^ 77 -r 



n=l 



OJV+I 



flAT + l 



E 



„ri-Af-l 
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Since all a„ and 6„ are positive, this is a decreasing function for x — >■ +00, and the 
lowest value of cLn will be the unique solution of 



M'-ES^ 



I \ -1 



N 



1 \{N + l)x^+^ 



n X" 



1. 



a.N+1 



We obtain the following table of values 



di 


1.03922 


de 


4.54145 


dii 


5.73661 


^20 


8.73298 


d2 


2.38568 


dr 


4.75734 


di2 


6.03061 


dso 


12.37349 


ds 


3.33232 


ds 


4.97336 


di3 


6.33969 


^40 


16.05983 


di 


3.92171 


dg 


5.20626 


du 


6.66091 


^50 


19.75448 


d5 


4.28707 


dio 


5.46090 


di5 


6.99175 


deo 


23.45053 



D 



Remark 5.6. The numbers dj^; in Lemma [53] are very similar to the numbers cjv of 
Lemma [5.31 This is no more than an experimental observation, but since the c^ 
numbers are easy to compute and d^ are somewhat elusive, it has been useful to 
start from cn as an approximation to d^ in order to compute d^- 

5.2. Some inequalities. 

Lemma 5.7. For u > 2 we have logali(u) < 21ogM. For u > e^ we have ali(u) < 
2m log w. 

Proof. The first inequality is equivalent to ali(M) < u^. Since li(a;) is strictly in- 
creasing, the inequality is equivalent to u < li(u^). 
For u > 2 we have h(u) > h(2) = 1.04516 ... so that 



\i{u^) = li{u) + 



"" dt 



\ogt 



>1 



logM^ 



Hence, our inequality follows from j" J^ > m — 1, that is from u > 21ogM. However, 
this last inequality is certainly true for u > 2. 

The second inequality is equivalent to u < li(2ulogu) and has a similar easy 
proof. D 



Lemma 5.8. For all integers n > 1 we have 
(48) 



e/- l0g"+4 

where /„ — 4(n + l)/3. 



" (log log i)" (log log u)"- 

dt < iu- 



log"+^ u 



{u> e-^") 



Proof. Notice that /„ > 1. For t > e the function log log t is positive and increasing 
so that 



r (log log i)" , „ , ,„ r dt 

/ \ „+i, dt < log log u)" / -^^ 

Jef^ l0g"+'i Je/„ l0g"+' 



It remains to be shown that 



dt 



< 



4u 



b„n+l I — 1 n+1 ' 
R t log U 



:g"-^ t 



Replacing u by e^ this is equivalent to 



^ e* , 4e^ 
dt< 



i"+i 



„n+l 



(u > e^"). 



(X > In)- 
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For the function 

n('r\ ■= — / 



4e^ 
G(x) ■= —— - / — — dt 






we have 

4(n+ l) 
x"-^^ \ X 

so that for a; > 4(n+ l)/3 we obtain G {x) > 0. Since G'(/„) > we have 0(0;) > 
for aU X > fn- □ 

Theorem 5.9. The polynomials Pn{y) defined in (j2ip satisy the inequalities 

(49) |P„(2;)|<3-n!y", y > 2, n>l 

and \Pa{y)\ <yfory>2. 

Proof. Since Po{y) = y — 1, the second assertion is trivial. 

Given r > 0, for each polynomial P{x) = X!n=o '^n^" ^^ define 



TV 
n=0 



1^11 = Ei«" 



It is easy to show that 

||P + QI1<I1P|1 + I1QI1, I1PQI1<I1P|1-|1QI1 

and that for the derivative of a polynomial of degree < N 

^ AT ^ N 

IIP'II = E n|a„|r"-i < - ^ |a„|r" ^ -||P||. 

n=0 n=0 

For y > r we have the inequality 

Af Af AT 

\P{y)\ = |5] a„y"| < J^ \a,,\y" < J] |a„|r"(y/r)" < {y/rf\\P\\. 

n— n— n— 

Hence, our Theorem follows if it can be shown that for n > 1 we have \\Pn\\ < 
3-2"n! (for r == 2). 

Define 5„ := ||-Pn||- By ()30p we have S*,! = — P„(— 2), and it can be shown that 
S™ <3-2"n! for <n < 15. 

For n > 15 it follows from (|26l) and the aforementioned properties of ||P|| that 

1 """'^ h 

Sn < nSn-1 + -zSn-1 + - ^ fcUfc - l)^*;-! + 5fe + - Sk-l) Sn-k-1 ■ 

k=l 

It follows that Sn < r„ where T„ ~ S„ < 3 • 3" n! for < n < 15, and that for 
n > 15 

3n Iv^V,^. K^k~2) \ 

Tn '■— -^Tn-l H 2^\ ^ ^ 9 Tk-ljTn-k-l- 

fc=l 
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Now we proceed by induction. For n > 15 and assuming that we have proved 
Tfc < 3 • 2*^ fc! for k < n, we obtain 







„ n— 1 

+ -E(fc2^ 


k\ + 


fc(3fc- 
2 


2)2fe- 


-i(fc-l) 


1^2"- 


-»-'(»- 


/c- 


1)! 


Hence 






















Tn 

3-2"n! 


- 4 


3 v-^ / k ■ k\(n 


-k 
■nl 


-'^ 


{3k 


-2)fc!(n 
8-n! 


-k- 


"^). 










3 3 

- 4 ^ 8n2 


>; 
fe=i 


7fc-2 


- 4 


3(7n- 
8n2 


»' + 


o "-2 

3 V . 

8n2 ^ 
fc=i 


7fc- 


-2 



Therefore, by using the symmetry of the combinatorial numbers, we obtain 



n-2 



Tn 3 3(7n-9) 3 ^7n-ll 



<- + i^—^ + -—Y,——i< 



3 • 2" n! - 4 8n2 16^2 ^-^ f""!) 

fc=l \ k J 



3 3(7n - 9) 3 ^, 7n-ll 

- 4 + 8n2 + 16^ ■ ^" " ^^ n-1 " 
3 3(7n-9) 3 ,^ , , n(4n - 63) + 87 , 

for n > 15. D 

Corollary 5.10. We have 

(50) |P„_i(j/)|<n!y", n > 1, y > 2. 

Proof. This foUows easily from the above Theorem. D 

5.3. Main inequalities. To simplify our formulae we introduce some notation. 
First we set r„ := 3 • nl so that, for n > 1, we have |P„(i/)| < r„y" when y > 2. 

Let c„ and d„ be the constants introduced in Lemma 15.31 and Proposition 15.51 
Let a„ be equal to max(e, c„, d„) and let /?« > e be the solution of the equation 

(51) = an- 

log a; 

(The function j-^ is increasing for t > e). 

Finally, define a;„ :— max(/3„, /„, e2), where /„ is defined in Lemma 15.81 

Proposition 5.11. Let x be a real number such that x > Xn, and set y :— logx. 
Then 

y > 2, a; > c„y, x > (i„y, x > fn- 

Proof. Since x > Xn — max(/3„, /„, e'^) we have x > e^ , so that y = logx > 2. 
We also have a; > /?« > e. Since j^^ is an increasing function for t > e we 

obtain j^ > jj|^ = a„ = max(e,c„,d„). Therefore, | > c„ and | > d„ as 
required. D 
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We insert a table of the constants 



Xi 


7.38906 


Xq 


10.81135 


Xll 


16.00000 


a;20 


29.57923 


X2 


7.38906 


X7 


11.70187 


2^12 


17.33333 


2^30 


47.86556 


X3 


7.38906 


Xs 


12.60164 


a;i3 


18.66667 


2^40 


67.69154 


X4 


8.29874 


Xg 


13.58167 


a;i4 


20.00000 


2^50 


88.57644 


X5 


9.77283 


2^10 


14.66667 


2^15 


21.42740 


2^60 


110.29065 



For each natural number N we set 



N 



(52) 



Wn = 1 + J2 



Pn-i{y) 



n=l 



and frequently we write W :— Wn when N is fixed. 

Proposition 5.12. For N > 1 let W = Wn (as in ([5^ ). Then for x > xn and 
y = log a; there exists with \0\ < 1 such that 



(53) 



W + xW + xWx + Wy - X - y - logW = e ■ rN+i 



y 



rN- 



Proof. Denote by T = T{x, y) the value of M^ + xW + xW^ + Wy - x - y -\ogW . 
Then we have 

T = (1 + a;) V '^"~^'^^'' - V '^^""^'^^^ + V I^niM- 

n— n— 1 n— 1 



From ([24)) we have the expansion 

(54) >o8(i+f:^)"^-E%l 



(») 



n=l 



Prom Proposition 15.111 we know that y = log a; > 2 and x > yd^- From (|50|) . for 
y > 2, we have \Pn-i{y)\ < nly" so that we have the majorant 



(55) 



N 



Pn-iiy) 



N 

<log(l-^ 



{ (x/yY 



(by considering this expression as a power series in a; ^, and y as a parameter). 
From dSll and (l55l). we obtain 



N 



(56) 



log 1 + E 



^n(y) 



N 

E 



Qn(2/) 



<S'Ar(x,2/) 



71—1 n—1 

where S]\[{^,y) is a power series majorized by the Taylor expansion of 



Mi-E 

(compare equation (j47| ). 



1 _ Y^ Qn 1 

^1 (a;/2/)"y ^ n (a;/?/)^ 
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By applying Proposition 15.51 we deduce that, for x > ydN, there exists 9 with 
\0\ < 1 and 

(57) ^-(-'^)-^^!^- 

If we substitute ([56|) in the expression for T, then by Theorem 14.31 ^H the terms 
in x^" with n < N cancel out, and the terms in a;~^ add up to — P/v(?/)x^^. It 
follows that 

(58) T = -^^ + SN{x,y). 

x" 

Therefore, since y > 2, we have 






N 



X 

so that from (|43)) , 

x™ V X / x" x" 

where TV > 1 and § + ^^^^ < 3 for a; > e^ are applied. D 



Proposition 5.13. For each natural number N let un — e^" . Then there exists 
vn > UN such that 

(59) HfNiu)) -u = 0- 13{N + 1)! ^(|Qg^Qg")'^ ^ (^ > „^) 

log u 



where \9\ < 1 and 



N 

(60) /Ar(e^) := a;e^W^Ar(a;,logx) = xe^(l + ^ 



F„_i(logj:) 



n=l 



Proof. To simplify the notation, the abbreviation W{x,y) = Wjv(a;,j/) is used. 
Differentiating (|60p we obtain 



^(li(/^(e^))-e-) 



^-——-(e'^W + xe'^W + xe"" (w^ + -Wy) ] - e^. 
Mle^)) I \ X J } 



l0g(/Ar(e-)) 



Assume that x > xat , so that x > d^ log x and x > e^ . We may apply ((53)) to 
obtain 



-(li(Me^)) - e^) = ^^^^^^{ly + .ly + .14^. + M.,} 



This may be simplified to 



'^ hUf (.^W -^-^ - ^"^ ^ T-jv+i log^a; 
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Since x > xn we have x > ycN, so that by Lemma [5?3l 

P„_.i(logx)^ / — "' 



<i + E^^^-f^)|^K^-E 



^1 i'^/yy' 



> 1 



that is xWn{x, \ogx) > 1, so that log(/Ar(e'^)) > x. Hence, for x > xn (with 
another 0), we have 

— (h(/^(e-)) - e) = 0^ e^ 

Defining Hn{u) :— li{fN(u)) — u the above equation is equivalent to 

, rN+i log^a: 
■^^(^ ) = ^ ^JvTi ' (^ - ^^) 

and, since um '■— e^™, 

WW ^ /j'^JV+i (loglogu)^ 

log u 

Lemma 15.81 can be applied since xm > Jn: so that u > un > e-^". Hence, 
integrating over the interval (un,u) we get 

Tj f \ u ( \ . a 4rjv+iM(loglogM)^ 

log u 
The function (A^+1)!t — tJ^^— -(loglogw)^ is increasing (as product of two positive 

increasing functions) for u > e^" , so that there exists v^ > un for which this 
function is greater than Hn{un), so that 

^ 13-(JV + l)!7.(loglogu)^ 
HnW^B , {u>vn). 

log M 

D 

Remark 5.14. For the values of n appearing in our tables, the equality Un — Vn 
holds, since, in these cases, 

„/ (ri + l)!u„(loglogu„)" 

l0g"+' Un 

Lemma 5.15. For any natural number N, and u > e^" we have log fN{u) < 
2\ogu. 

Proof. First observe that the hypothesis u > e^" implies (with u — e^) that x > 
Xn, so that logx > 2 and x > cjm logx. fProposition lS.lip . 

The inequality log /at (u) < 21ogu is equivalent to fN{u) < u'^, and together 
with u = e^ it is equivalent to 

Pn~lil0gx)\ ^ 2:e 



-(l+E 



xe-^ I 1 + > ^^^^^^ I < e^ 



From Corollarv l5.10[ since x > e^ and x > cat?/, and by Remark l5.14[ 
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Hence our inequality follows from 

X N I 

ye X ( v-^ n\ \ 

Since we assume that x > ycN, the second factor is greater than 1, so that 

X y V _ {x/y)^ J X 

Finally, it is easy to prove that e^ log a; > x^ for x > e^. D 

The asymptotic expansion with bounds can now be proved. 
Theorem 5.16. For each integer N > 1 

(61) ^l[(u)^fj,{u) + 260iN + iy.u(^-^f^^^^Y, {u>vm), 

\ log u / 

where vn is the number defined in Proposition \5.13\ . 

Proof. Since li(ali(u)) = u, Proposition 15 . 1 3l vields. for u > vn, 

liiMu)) - h(an(.)) = / -£- = miN + 1)! u^ip^l?^. 

Jali(«) fog I log ^ U 

Since vn > un = e^", u > vn implies logw > 2, hence u>2. 

From Lemma 15.71 logali(w) < 21ogw, for u > 2. Analogously, Lemma 15.151 
implies that log /at (m) < 21ogu, for u > e^™. Therefore, for u > ujv, we have 



|ali(u)-/Ar(u)| ^ 



/iv(u) ^^ 
ali(«) log* 



21ogu 

It follows that there exists 9' with \9'\ < 1 such that 

r/iv(u) ^^ 



ali(«)-/jv(u) = 0'(21ogu) 

Ja,li(u) ^Ogt 

and the result follows easily. D 

The actual error appears to be much smaller than that given in Theorem 15.161 
However, as usual with asymptotic expansions, having a true bound allows realistic 
bounds to be given of the remainder for specific values of N. 

The true error after N terms of an asymptotic expansion, while the terms are 
decreasing in magnitude, is often of the size of the first omitted term. In our case, 
the magnitude of the term P/v(logx)a;~^~^ depends on the polynomial PN{logx). 

Numerically, it appears that for n > 3: 

(62) \Pn{y)\<(-, ) y", (y>21ogn) 

V e log n / 

although we have not been able to prove this. 

From Theorem 15.161 more realistic bounds can be obtained for the first values 
of N. This is done in the following Theorem. 
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Theorem 5.17. For 2 < N <11, we have 

n=l 



(^^) ^^ = 1 + 2^ -n +^-20-i^to^J •^^^' ^">"^) 

where 



Z2 = 1.50, Z3 = 2.34, Z4 = 3.32, zg = 4.33, zg = 5.36, 
Z7 = 6.39, zg = 7.43, zg = 8.46, zio = 9.50, zn = 10.53. 

Proof. By taking A^ = 10 in Theorein l5.16[ we have, for u — e^ > e^'-", (recall also 
Remark ??) 

ali(e^) -^ P„_i(loga:) log^°x 

^^ = 1 + 2^ ^ + ^^^T^' (^>2:io) 

n—1 

with i? = 26 • 11! = 1 037 836 800. 

We compute the maximurr|3 A/„ of |P„_i(loga;)/log"^^ a;| for a; > xio, so that 
for any 2 < A^ < 10, we have 

ali(e"^) ^ -^ ^ y^ P„-i(loga:) ^ 



log^/ f. P„_i(log:r)log"-^^-^x 1 



n=Af+l 



SO that 




aH(e^) 


^ 1 y> P„-i(loga;) ^ 




, ,log~x^ ^ Af„log"-^-^x ^ „W"-^x 




n=Ar+i 



We determine a value z^ > xio such that, for x ~ z^, 

n=Ar+l ° 

Since this is a decreasing function of a;, we obtain for x > z^ 

(64) «:l = l + V^i.iMfi)+..20(-^)"i^. 

n— 1 ^ 

We consider the function 

N 



ali(e-) ^_y.P„-i(logx)^ a;^+i 



on the interval (1.3, z^), to determine the least value of zm for which (|M|) is true. 



3a/2 = 1, A//3 = 1/2, M4 = 1/3, Afs = 0.250636, Afg = 0.526887, Afy = 1.300565, Afg = 
3.719653, Afg = 12.070813, Afio = 43.788782. This last maximum would be much smaller if the 
maximum were taken from a point slighthly greater than xiq. 
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In this way we find: Z2 = 32 and then Z2 = 1.5; Z3 = 49.5 and then Z3 = 2.3395; 
Z4 = 82 and then Z4 — 3.3114; Zg — 155 and then Z5 = 4.3237. 

If we take TV = 20 in Theorem WM. we obtain z^ = 113, z'-j = 143, Zg = 187, 
z^ = 251 z[q = 353, z^i = 528 from which zg = 5.3514, Z7 = 6.3851, zs = 7.4208, 
zg = 8.4566, zio = 9.4914 and zn = 10.5251 are obtained. D 

Remark 5.18. We have proved ((63|) only for 2 < iV < 11, ahhough something 
similar appears to be true for the general case. If (|63p were true for all n, then for a 
large u = e^ we could take N fn x terms in the expansion and in this way the error 
would be « -2%, so that ali(u) could be computed with an error less than sa 20. 

In fact, for several values of u, the terms of the expansion have been computed 
up to the point where these terms start to increase. Always the computation is 
terminated when N ^ x and the error appears to be bounded. ( For example, 
with u — 10^'^*', we compute 230 terms of the expansion, which coincides with 
loglO^'^'^ « 230.259. The approximate value obtained for ali(u) has an absolute 
error equal to 40.94738, which can be compared with the fact that ali(10^°°) has 
103 digits ). 

6. Applications to p„. 

6.1. Asymptotic expansion of p„. Inequalities for the n-th prime number can 
be found in [18], [16], [12], [4]. In fact, from 7r(a:) = li(a;) + 0{r{x)) we may obtain 
Pn = ali(n) + 0(r(n log n) log n), if r{x)/x is sufficiently small. For example, in 
[l2] , it is noticed that from a result of Massias [10] , it follows that 

(65) Pn = ah(n) + Oine-"^^"^) 

so that the asymptotic expansion of ali(n) is also an asymptotic expansion for p„, 
that is, 



N 



Pfc_i(loglogn)\ / /log log nN^ 



V ^ log^n ) V V logn 

By assuming the Riemann hypothesis, Schoenfeld [H] has proved 

(67) |7r(a;) - X\{x)\ < — V^logx, {x > 2657). 

87r 

This result will be used to obtain (under RH) some precise bounds for p„ 

Lemma 6.1. We have -^^(loga::)^ < ali(x) for x > 94. 

Proof. The inequality is equivalent to 

h(^/a;(loga;)5) < x. 

Differentiating the function f{x) :— x — li(-/i(loga;)5) we get 

^v^^_l 1 ^ log'^'^ , 51og^/^ x^ 

log(Vi(logx)i)V 2^ 2Vi 

Hence this derivative is positive if and only if 

, 3/2 /logx 5\ ^/loga; 5, , 
log3/2a;|^_|_ + _j <7^(^_|_ + _ioglogx 

For a; > 94 we have (log x) 2 < y/x, so that /' {x) > for x > 94. 

Finally, one may verify that /(94) > 0. Hence f{x)>0 for a; > 94. D 
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Theorem 6.2. The Riemann hypothesis is equivalent to the assertion 

1 5 

(68) \pn — ali(n)| < —\/nlog'^ n for all n > 11. 

TT 

Proof. First we assume the Riemann Hypothesis and prove (|68|) . Let r(x) := 
g^^Y^loga;, f{x) := \i{x) — r{x), and g{x) := h(a;) + r{x). For x > 1 we have 
f{x) < \i{x) < g{x), where the three functions are strictly increasing. From (j67p 
for X > 2657, we also have f{x) < tt{x) < g{x). 

The inverse functions satisfy g^^{y) < ali(y) < f^^{y), and if y = n > 7r(2657) — 
384 is a natural number, then g~^{n) < p„ < f^^{n). It follows that the distance 
from ali(n) to p„ is bounded by 

\pn — ali(n)| < max(/~"^(n) — ali(n), ali(n) — g~^{n)). 

Hence, we have to bound f~^{y) — ah(2/) and ali(y) — g~^{y). 
We consider y as a parameter and set a — ali(y), so that li(Q;) — \\{a\\{y)) — y. 
Consider the function m(^) := /(^) — li(a) = f{0^ Vi which is strictly increasing 
and satisfies 

uiO=im-r{0-Ho.)= (^ ^-r{i). 



Therefore, u[a) — —r{a) < and 

u{f-\y)) - f{f-\y)) - li(a) = y - h(ah(y)) = 0. 

If apoint b is found where u{b) > 0, then a < f^^{y) < b, so that b—a > f^^{y)~a 
and one of the required bounds is obtained. 

Therefore, we try b = a + c^{\ogy)2 with c < 1. We have 

u{a + Cy/y {log y)^ ) ^ / j — - - r{a + Cy/y{\ogy)^) > 

> , , ^7^ ;t- - r{a + cV^(log y)^). 

log{a + Cy/y{\og y) 2 ) 

From Lemma |6T] for y > 94, we have y/y{logy)^ < ali(j/) — a, so that 



^a + cVy(logy)i)>^^f|^-K2a) 



c^i\ogy)-2 ~ ^V2^\og\2a) 



log(2a) 

We want to show that this expression is positive. For y > 94, we have a — ali(y) < 
2y\ogy (by Lemma [??7|) . so that a < 2y\ogy < Ay logy < y^ (for y > 94), which 
yields (with c = l/vr) 

5 I , — 2 
Cy/y{logy)^ - — V2alog (2a) > 

OTT 

> Cy/y {log y)^ - — y^4y log y log^ (4?/ log y) > 
ovr 

> cVy(log2/)^ - — V4ylogylog^(2/^)) = 0. 

ovr 

Hence, we have proved that f^^{y) — a < ^ ^yy (log yY^^'^ for y > 94. 
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To bound a — g^^{y), we consider the function v{^) :— g(^) — li(a) = 5(0 ^ V- 
Then 

and v{a) ~ r{a) > 0, v{g^^{y)) — g{g^^{y)) — y — 0. If a value b is found such 
that v{b) < 0, it wih follow that a — g~^{y) < a — b. 

Choose 6 = a — Cy^(logy)^/^ with c = i. We claim that v{b) < 0. We have 

v{b)=r{b)- r-^<r{b)~^ 
Jb logi log a 

and our claim will follow from r(b) log a — c-^y (log y)^^^ < 0. Finally, since b < a = 
ali(y) < 2ylogy, and by Lemma ISTZl 

1 r^T^ o 1 



r(5)loga < —y/ali{y)(logal[{y)f < —y/2y\ogy{2\ogyf 

OTT OTT 

which proves our claim. 

Hence, (by assuming RH), we have proved that |p„ — ali(n)| < i-y/n(logn)2 for 
n > 384 > 94. By verifying all 1 < ti < 385, we find that the inequality holds 
except for n < 11. 

The reverse implication is simple. From |p„ — ali(ri)| = 0{n^^^) for any e > 0, 
we may derive that 7r(x) = li(x) + 0{x^^^). It is well known that this is equivalent 
to the Riemann Hypothesis. D 



Remark 6.3. The inequality (j68p is only proved by assuming the Riemann Hypoth- 
esis, but is stronger than those contained in [18], [16], [12], [4]. Inequality (|68l) gives 
approximately half of the digits of p„. If our conjecture that ([63]) is true for all N 
is also assumed, then the asymptotic expansion gives about half of the digits of p„. 

6.2. Inequalities for the n-th prime. Let 

f \ 1 /i , Y^ Pfe-i(loglogn)\ 
sn[ti) = nlogn 1 + > — v 

where sq — nlogn. 

Cipolla noted that for fc > 1, Pu{y) = (-l)''+^%- + - • • , andPo(2;) =J/-1- Hence, 
except for the first term, eventually the sign of the fc-th term Pfe_i (log log n) log^ n 
becomes (—1)'^. The asymptotic expansion implies that there exist rj.^ such that 

p„ > so(n), n > ro, p„ > si(n), n > ri, 

Pn<S2N[n), n>r2N, Pn> S2N+l[n), n>r2N+l- 

In fact, ro = 2 is the main result in [IS], ri = 2 is proved in [4] and r2 = 688 383 is 
proved in [5]. The value of r^ ioi N > 3 has not been determined. See Theorem 
16.41 for an estimation of r^ by assuming RH. 

The above reasoning may give the impression that the terms of the asymptotic 
expansion of ali(w,) are alternating in sign, starting from the second term. However 
this is not true. For example, computing the first 230 terms for ali(10^°°), we found 
only three positive terms Pq/x, Pi/x^, and P^/x'^. In fact, the sign of the fc-th term 
is that of Pfc-i(loglogn). Thus we are interested in the sign of these polynomials. 
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The polynomials P/v(j/), for 1 < iV < 23 of odd index, have one and only one 
real root, which is positive. Starting from Pi{y) which vanishes at j/ = 2, these 
roots are 

2, 4.23415, 5.83131, 7.43591, 9.07979, 10.6881, 12.2538, 
13.7876, 15.2977, 16.79, 18.2683, 19.7353. 

The polynomials P2, P4 and Pq have no real roots, and all Pg, . . . , P22 have two 
positive real roots. These pairs of roots are: 

(6.4306,8.2185), (7.16158,9.88528), (7.90293,11.4752), 

(8.63359,13.0241), (9.3507,14.5452), (10.055,16.0458), 

(10.7478, 17.5307), (11.4307, 19.003). 

For example, P9(loglogn) is positive only for n > exp(e^''^'), which is a very 
big number. 

The even terms at first sight appear negative. However Pio(loglogn), for exam- 
ple, is negative except in the interval exp(e^'^^') < n < exp(e^'^*'). 

In a certain sense, the inequalities (|69| are the wrong inequalities. These in- 
equalities would hold only for very large values of r^, especially when we want a 
lower bound of p„ > S2N+1 (except for the three known cases). We estimate r^. 

Theorem 6.4. Let r^ be the smallest number such that 

log log n — 2 



(70) p„ > S3 := 7ilogn + ?i(loglogri - 1 



logn 

(loglogn)^ — 61oglogn + 11 
21og^n ' 



n > r^. 



Then, if the Riemann Hypothesis is assumed, 

39 X 10^9 < ra < 39.58 x 10^^ 
Proof. By Theorem 16.21 there exists 9i with 16*11 < 1 such that 

p„ = ali(n) -1- 6*1 -^- — (log 71)2, n > 11. 

TT 

By Theorem I5.17[ with 5 < iV < 10 and setting n — e^, x — logn, and y 
log log n, we have for n > e^™ 

y-1 y-2 y2_6y + ll 



aU(n)=a;e^fl + ^^ + 

V X 



X X 2x^ 



N 



Psiy) , ^Pfc_i(y) y 

1 h f 2CJV 



E 



N 



fc=5 

The inequality of the Theorem is obtained if 

V X^ ^-^ X'^ X"^'-/ TT 

k=5 

This is equivalent to 

P3(2/) 4-5:Pfc-i(2/)e-('=-4)^ > c^y^e-(^-3)^ + ie^^e"^^". 

k=5 
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Since ^3(2/) — > +00 and all the other terms tend to as y — > +cxi, it is clear that 
the inequality is true for y > yo. With N = 10, we find yo = 4.254946453 . . . The 
inequality p„ < s„ is true for n > 3.95702224148845656 x 10^°. This proves that 
rg < 39.58 X lO^^. 

In order to show that rg > 39 x 10^^, we directly show that, for n = 39 x 10^^, 
the opposite inequality p„ < S3 is obtained. 

We compute 

S3 = 2.87527 18639 02974 79681 42399 35057 89294 02005 87915 X lO^^ 

Now we can compute ali(n), for which we already have obtained a good approxi- 
mation through the asymptotic expansion, and then apply the Newton method 

ali(n) =2.87527 18639 02495 21516 14800 14732 45414 39731 X 10^^ 

Therefore, from Theorem 16. 2[ we obtain p„ < ali(n) + i^/nlog^ n, so that 

Pn < 2.87527 18639 02756 97808 39055 05640 30082 86370 11482 X 10^^ 

and we can conclude that p„ < S3. D 
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1. Introduction. 



C " 3 Abstract. A new derivation of the classic asymptotic expansion of the n-th 

^Nj , prime is presented. A fast algorithm for the computation of its terms is also 

given, which will be an improvement of that by Salvy (1994). 

O , ' Realistic bounds for the error with li~^ (n), after having retained the first m 

J^ ^ terms, for 1 < m < 11, are given. Finally, assuming the Riemann Hypothesis, 

^^ , we give estimations of the best possible ra such that, for n > rs,, we have 

Xf^ ' Pn > S3(n) where s^in) is the sum of the first four terms of the asymptotic 

expansion. 

H 

C^ I l-l- Historical note. Chebyshev failed to fully prove the Prime Number Theorem 

(PNT), but he obtained some notable approximations. For example, he proved that 
for every natural number n: if the limit 

<y' lim i^i^(4x)-li(x)) 

~^ , exists, then this limit must be equal to 0. 

.^ ■ The question was decided by de la Vallee Poussin (1899) when he gave his bound 

\f^ , on the error in the PNT: The above limits exist and equal 0. 

In 1894, Pervushin, a priest in Perm, published several formulae obtained em- 
("^ , pirically about prime numberqj. One of these formulae gives the following approx- 

^SJ ' imation to the n-th prime 



— = log n + log log n - 1 + — — 5— . 

n '^ '^ ^ 121ogn 241og^n 

Cesaro then published a note [1, 1894] where he asserts that the true formula is 



H 

•^ ■ P« 1 , 1 1 1 , loglogn-2 

■ - - — = log n + log log n — 1 



n log n 

(loglogn)2-61oglogn + ll , ,, _2 x 

— ^ + o(log n). 

2 log n 

Despite no mention by Cesaro in [1], the editors of his collected works added a 
note to \Xi pointing out that certain formulae quoted by Cesaro, since they followed 
from the results of Chebyshev, were only established under the assumption of the 
existence of the implied limits. It therefore remains unsurprising that Hilbert, in 
the Jahrbucto stated that Cesaro did not prove his formula. 



Ivan Mikheevich Pervushin (1827-1900) (HBan Mhxccbhh IIcpByniHH). No small achieve- 
ment if we note that he had only a table of primes up to 3 000 000. 

Jahrbuch iiber die Fortschritte der Mathematik (1868-1942), a forerunner for the Zentralblatt 
fiir Mathematik, at present digitalized at http://www.emis.de/MATH/JFM/JFM.html 

1 
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Landau J7, 1907] several years later was better informed: a formula, like that of 
Cesaro, would imply the PNT, which had yet to be proved at Cesaro's time. How- 
ever, using the results of Chebyshev, Cesaro may claim that if there is some formula 
for Pn correct to the order n(logn)~^, then it must coincide with his formula. 

Cipolla 1902] obtained an infinite asymptotic expansion for p„ and gave a 

recursive formula to compute its terms. He published after the results of de la 

Vallee Poussin but it seems that he was unaware of these results, so that gave 

his proof under the same hypotheses as Cesaro. So uninformed was he that he 

attempted to prove some false formulae of Pervushin already corrected by Torelli 

[22] 

loglogn-1 ,^/loglogn\2 

Pn+i - Pn = log n + log log n H h e> — 

logn V logn / 

with an impeccable proof that if such a formula exists, then it must be this formula. 
(Such a formula would refute the twin prime conjecture, and today the above 
formula is known to be false.) 

In his Handbuch Landau [SI § 57] obtained by means of the procedure of Cesaro, 
some approximative formulae for p„, and explained that the method could give 
further terms. He also mentioned some recursive formulae without giving any clue 
for their derivation. 

We may say that Pervushin was the first to deal with a formula for p„, albeit 
that he gave only the first few terms. Cesaro then proved that in the case such 
a formula exists, it must be one from which he would be able to derive several 
terms. Cipolla found a method to write all the terms of the expansion if there is 
one. Landau saw that the results of de la Vallee Poussin imply that the expansion 
certainly exists. 

The algorithm given by Cipolla is not very convenient for the computation of the 
terms of the expansion. He iteratively computes the derivative of some polynomials 
appearing in the expansion but computes the constant terms as determinants of 
increasing order. Robin [ITl 1988] considers the problem of computing these and 
other similar expansions, leaving the problem of computing the constant terms of 
the polynomials as an open problem. Later Salvy [20l 1994] gives a satisfactory 
algorithm. This algorithm needs O{n7'^y/\ogn) coefficient operations to compute 
all the polynomials up to the n-th polynomial. 

The asymptotic expansion of p„ also plays a role in the study of g{n), which is 
the maximum order of any element in the symmetric group 5„. In fact, log g(n) 



has the same asymptotic expansion as y li (x) |llj . 

There are many results giving true bounds on p„, for example we mention p„ > 
nlogn [18J 1939], and p„ > n(logn + log logn - 1) [41 1999] both for n > 2 (with 
partial results given in [19], [16], [12], [4]). In [5] it is also proved that 

/ log log n — 2 \ 

Pn < n\ log n + log log n — 1 H , n> 688 383. 

V log n / 



1.2. Organization of the paper. In this paper we present a new derivation of 
the asymptotic expansion for p„ and obtain explicit bounds for the error. 

First, it must be said that the asymptotic expansion has, in a certain sense, 
nothing to do with prime numbers: it is an asymptotic expansion of ali(x) := li~ x 
which is the inverse of the usual logarithmic integral function. 
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In Section [3] a proof of the existence of the expansion is given, fohowing the path 
of Cesaro, since it cannot be found elsewhere, ahhough it is frequently claimed it 
can be done. This Section is not needed in the rest of the paper. 

In Section|4l a new formal derivation of the expansion is given. We obtain a new 
algorithm to compute the polynomials (Theorem 14.91) . This is simpler than that 
given by Salvy [20' . Our algorithm allows all the polynomials up to the n-th one to 
be computed in O(n^) coefficient operations (Theorem 14. 11 1) . It must be said that 
these polynomials have 0{'n?) coefficients. 

In Section [SJ independently of Section [31 we prove that the formal expansion 
given in Section |4] is in fact the asymptotic expansion of ali(a;) and gives realistic 
bounds on the error (Theorem 15. 161 and 15. 17|) . 

In Section [6l the results are applied to p„ the n-th prime. Using the results of 
de la Vallee Poussin it can be shown that the asymptotic expansion of ali(ri) is also 
an asymptotic expansion for p„ . 

By assuming the Riemann Hypothesis, we found (Theorem I6.2p that 

\pn- &\\{n)\< -y/n{\ogn)i , n > 11. 

TT 

This bound of p„ is better than all the bounds cited above. 

We end the paper by motivating why the above bounds have not been extended 
to further terms of the asymptotic expansion (Theorem 16. 4|) . 

Notations: With a certain hesitation we have introduced the notation ali(a;) to 
denote the inverse function of li(a;). 

In Section [51 where explicit bounds are sought, it has been useful to denote by 
9 a real or complex number of absolute value \9\ < 1, which will not always be the 
same, and depends on all parameters or variables in the corresponding equation. 

Acknowledgement: The authors would like to thank Jan van de Lune ( 
Hallum, The Netherlands ) for his linguistic assistance in preparing the paper, and 
his interest in our results. 

2. The inverse function of the logarithmic integral. 
Usually li(a;) is defined for real x as the principal value of the integral 

h(x) = p. V. / . 

^ ' Jo log^ 

It may be extended to an analytic function over the region Q ~ C\ (— oo, 1], which 
is the complex plane with a cut along the real axis x < 1. The main branch of 
the logarithm is defined in fl and does not vanish there. Therefore, li(z) may be 
defined in fl by 

(1) ii(^)^ii(2)+ r^, zen 

where we integrate, for example, along the segment from 2 to z. 

For real a; > 1, the function li(a;) is increasing and maps the interval (1,+cx)) 
onto (— oo, +oo), so that we may define the inverse function ali: R — )■ (1, +cx)) by 

(2) h(ah(x)) = X. 
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The function li(2:) is analytic on f2, so that ali{x) is real analytic. It is clear that 
we have the following rules of differentiation 

(3) — \i{x) = , — ah(a;) = logah(a;). 

ax log X ax 

It is well known that the function \i{x) has an asymptotic expansion: 
Theorem 2.1. For each integer N > 

logxV ^log'^a; Mog^+ia;yy' 

This may be proved by repeated integration by parts (see [13l p. 190-192]). 

3. Asymptotic expansion of ali(a;). 
In this section, we prove the following 
Theorem 3.1. For each integer N > 

n—l 

where the P„_i(z) are polynomials of degree < n. 

In the case of A^ = the sum must be understood as equal to 0. 
The theorem says that, for each N, there exists an xn > 1 and a constant Cn 
such that 

ah(e^) ^ .^P„_i(logx) '-^+1 



V^ -r„-Uiogx-j log X 

1-Z^ Z^ ^^^ ^N+l ' {x>xn). 



xe-^ 

n— 1 

In the course of the proof we will make repeated use of the following 
Lemma 3.2. Let f{x) be a function defined on a neighbourhood of x = such that 
(6) f{x)^aix + --- + aNX^ + 0{x^+'^), {x -^ 0) 

where the ak are given constants. Assume that g{x) satisfies 

N n \ 1 N+1 

f7\ ( \ \^ Pni^Ogx) .^f^Og_X\ t ^, \ 

n—l 

where the Pn{z) are polynomials of degree < n. Then there exist polynomials qk{z) 
of degree < k such that 

^ n \ 1 N+1 

(8) f^g(x))^^^—^ + 0^ j, (X^+OO). 

Proof. It is clear that, for each 1 < n < A we have p„(loga:)a;~" = O((loga:/x)"). 
Therefore, g{x) = O{\ogx/x) (this is true even when A = and there is no p„). It 
follows that linij; _>._(_ 00 .9(2;) = and by substitution in ([5]), we obtain 

N -, Af+l 

(9) /(g(x)) = 5]a„g(x)" + 0^°^ "^ 
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By expanding the powers g{x)^ by ([7]) it is easy to obtain an expression of the form 

^ n \ 1 N+1 

(10) gix) ^^ -^ + e)(^___j, (x^+oo) 

fc=l 

where each Pn.k{z) is a polynomial of degree < k. By substituting these values in 
equation ^ and collecting terms with the same power of x, ([5]) is obtained. D 

We will prove Theorem 13.11 by induction. The following theorem yields the first 
step of this induction. 

Theorem 3.3. 

(11) ^^,^aC-^^), {X-.+00). 

a; log a; V logx / 

Proof. From Q with iV = we have "'^\'°g^ = 1 + O(log"^ y) for y ^ oo. Since 
limj;^_i_oo ali(a;) — +oo we may substitute y ~ ali(2:) and obtain 

, , X log ali(a;) 

(12) ° ^ =1 + 



ali(a;) Vlogali(.T) 

By taking logarithms 

log X - log ah(a;) + log log ali(x) = 0( -— — 

Vlogali(x) 

we obtain 

log a; _, , ,^/loglogah(x; 



logali(2;) V logah(a;) 



(13) r^ = i + K- 

logah(a;) V 

and it follows that 

(14) lim ^^ = 1. 

x-)-+oo log ali(x) 

By taking log in P^ 

'log log ah (a;) 



log log X — log log ali(a;) — 0(- 



log ali(x) 



we obtain 


loglogx I \ oi ^ 




loglogali(x) ' "'Vlogali(x) 


so that 




(15) 


log log a; 
Jl™ ^_,_^,^, , =1. 



In view of p4)) and p5)) . we may write (|12l) and (1131) in the form 

^^g^ xlogah(x) _^ ^ gf M, log^ _i I ^>;>^loglog^ 



ali(a;) Vloga;/ logah(a;) V log a; 

and by multiplying these two, we obtain 

xlogx _ _^ ^^loglogx 



ali(a:) V logo; 

from which (jlip can easily be deduced. D 
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Proof of Theorem \3.1\ We proceed by induction. For A^ = 0, our theorem is simply 
Theorem 13.31 with e^ instead of x. 

Hence we assume ([S]) and try to prove the case A'^ + 1. 

Our objective will be obtained by starting from the expansion of li(y). By ([4]) 

By substituting y = ali(e^) and applying (ITi)) we obtain 

(17) e-logali(e-)^ ^^ k\ ^o(^—\ 

^ ' ah(e^) ^ (logah(e^))fc Va;^+2y'- 

From our induction hypothesis, the expansion of logali(x) and (logali(a;))~'^ is 
now sought. 

By taking the log of ^ we obtain 

,og.l.(.') = . + log. + log{l + £ ^l^^fifi + 0(^1;;;^;)}. 

n— 1 



Lemma 13.21 may be applied with 



log(l + X) = X-^ + ^~--- (-1)^+% + 0(X~+i) 



-^ Q„+i(log.x) , ^^/log^+^x^ 



to obtain 

logali(e )=a; + logx + ^ — + 0[^ ^^^^ j. 

n—l 

The reason why we have written Qn+i instead of Qn{x) is revealed below. The 
above may be written as 

logah(e)^x|l + ^ + )^ ,;»+i +^(^lvT^J| 

Tl=l 



(18) logaH(e^) = a;{l+ ^ 



^+1 ^ n,.„^^ .i^„A^+2 



O"(loga;) , ,^/log™+^x 



O 



X" 
n—l 



rN + 2 



where the Q„(z) are polynomials of degree < n. Observe that knowing the ex- 
pansion of ali(e'^) up to (loga;/a;)''^+^ has enabled us to obtain logali(e'^) up to 
(logx/x)^"'"^; this will be of great importance in what follows. 
From (|18p . for all natural numbers n, 





1 1 


{:. 


N+1 

Tr 

fc=l 


9/c(l0g 


^'+c.(i 


og^+^. 




(logaH(e^))" x" 


^N+2 


By 


applying Lemma |3.2| with 




Af + 1 


/ Tl\ 








(1 + =^)-" 


- 1 = 


r=l 


I :) 


x' + 0{x 


-N-2\ 



)}' 
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we obtain 



1 1 / 

lloffalife^)!" ~ x" V ^ 



^1 T/ .n„„^^ ,i^„A^+2 



Ka(loga:) , ^(\og''+'x 



{logali(e^)}" x" V ^ x*^ V x^+'^ 

where the Vn.k{z) are polynomials of degree < k. 

By substituting these values of {logali(e^)}^" in ([TT]), we obtain 

e"logali(e") ^ , , ^' Ukjlogx) / log^+^x x 

ah(e^) ^^ x'= V x^+2 )' 

Hence again from ([T51) with n = 1 

-'V+l TA /I N 1 JV+2 

— T4,fe(log2;) „/log ^ a;^ 



ali(e^ 



i{,,j;Vi^,,(i^)J 



fc=l 



-'V+l TT n \ 1 -/V+2 

{i+i:^^+c>(i2|^)} 



fe=l 



from which we derive that there exist polynomials Wk{z) of degree < k such that 
xe- _ ^+^Wk{\ogx) , ^^/log^+'' 



^ ^ Wk{^ogx) .^/ iog x \ 

^+Z^ X'= V X^+2 y*- 



ah(e-) 



Another application of Lemma 13.21 yields 

ali(e-)^ ^^ T,(logx) / log^+^x 

xe^ ^ xfe V x^+2 

fc=i 

with polynomials Tk{z) of degree < k. Therefore, we have an asymptotic expansion 
of type ([Sj with A^ + 1 instead of N . The usual argument of uniqueness of the 
asymptotic expansion applies here so that Tk{z) = Pk{z) ior 1 <k < N . D 

4. Formal Asymptotic expansion. 

First we give some motivation. We have seen that the asymptotic expansion of 
ali(e^) is 

an(e^) = xe^F(x, log x), where F(x, y) := 1 + ^ ^n 

n—l 

and differentiation yields 

e^ logah(e^) = (e^ + xe'')V + xe^K + e'^Vy. 
Here logali(e^) — log(xe^y(x,logx)) = y + x + logV^, so that 

y + x + logV = V + xV + xVx + Vy 
which we write as 

(20) V ^l+'^--V-V^--Vy + -\ogV. 

XX XX 

This ends our motivation for considering this equation. 
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Consider now the ring A of the formal power series of the type 

y^ Qnjy) 

where the q-niy) are polynomials with complex coefficients of degree less than or 
equal to n. In particular go(y) is a constant. 

It is clear that A, with the obvious operations, is a ring. The elements with 
go = form a maximal ideal /. An element 1 + m with go = 1 is invertible, with 
inverse 1 — u + u^ — ---. It follows that if a ^ /, then a is also invertible. Hence / is 
the unique maximal ideal and A is a local ring. If a € A is a non- vanishing element, 
then there exists a least natural number n with g„(y) ^ 0. We define deg(a) = n 
in this case, with deg(O) — oo. 

As is usual in local rings, (see [9]) we may define a topology induced by the 
norm ||a|| = 2^'^'^^^°-\ which, with the associated metric, induces a complete metric 
space. Indeed A is isomorphic to C[[X, Y]\, by means of the application that sends 
X !->• x~^ , Y i-> yx~^. 

Given a G A with a — Yll^=o •^ , we define two derivates 



= 1 -^ n=l 



X" 



Finally the set U C A oi elements with qo = I form a multiplicative subgroup of 
A* (the group of invertible elements of A). For 1 + u E U, we define 

oo J^ 

k=l 

which is a series that is easily shown to converge since u'' E I^ . 
We are now ready to prove the following 

Theorem 4.1. The equation (|20p has one and only one solution in the ring A. 

Proof. For V^ G [/, we define T{V) as 

T{V) ■.= l + ^--V-V.,^-Vy + - logK 

XX XX 

It is clear that T{V) € U . We may apply Banach's fixed-point theorem. Indeed, we 
have deg(r(y)-T(VF))<l+deg(F-'l^), so that \\T{V)-T{W)\\ <\\\V-W\\. 
By Banach's theorem there is a unique solution to F == T{V). We may obtain 
this solution as the limit of the sequence T"(l). In fact, since deg(T(F) — T(W^)) < 
1 + deg(y — W), in each iteration we obtain one further term of the expansion. In 
this way, it is easy to prove that the solution is 

w- 1 w-2 

V = 1 + ^ + ^-5- + --- 

X X 

However, we are going to find more direct methods to compute the terms of the 
expansion. D 

Definition 4.2. Let V be the unique solution to equation (|20p . Since it is in A, it 
has the form 

(21) n-,j^)-i+f:^^ 

n— 1 
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where for n > 0, Pn-i(y) is a polynomial of degree < n. 

In the following sections, we prove that V yields the asymptotic expansion of 
ali(e^). For this proof the following property is crucial. 

Theorem 4.3. For N >1, let 

(22) w{x,y)=WN{x,y):=l + f2^^^^- 

k=l 

Then 

Proof. By the definition of the P„ we know that deg(T^ — W) > N + 1. Therefore, 
deg(T/ - T{W)) >N + 2. That is 



V-T{W) 



"0(2/) , vo{y) 



X 



N+2 ' rf.N+3 



where wq, vq are polynomials. We also have 



so that 



That is, 



V 


= w + 


00 

>: 


Pn 


-1(2/) 








uoiy) 

^N+2 + 


Voiy) 

xN+3 


+ ■ 


n 


00 

>; 

,=Af+l 


Pn-1 




X" 



W-T{W) 



w-Tiw)= ^^^^^ ■ "^^^ ■ ""^y^ 



j,Af+l ' x^+2 ' x^+3 

for certain polynomials u, v, . . . D 

In the sequel V will denote the unique solution to (PO)) . The element logF 
belongs to A, so that there are polynomials Qn{y) of degree less than or equal to 
n such that 



(24) fog^-E 



Qniy) 



X" 
n=l 



From equation (|20p . we may obtain logV^ in terms of V and its derivatives. It is 
easy to obtain from this expression the following relation 

(25) Q„(2;)=P„(2/)-(n-l)P„_i(y) + P,;_i(y), (n > 1). 

Theorem 4.4. The polynomials Pn{y) that appear in the unique solution (|2ip to 
equation (|20p may be computed by the following recurrence relations: 



Pq ~ y — 1, and for n > 1 

(26) ^ ""' 

i-„ = ni-n-i - n,-i + 

n 



-. n— 1 

p„ = 7ip„_i - p;,_i + - ^ fc{(fc - i)Pfc_i - Pfc - p^^_ jp„_fe_i. 



fc=i 
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Proof. By differentiating ([M]) with respect to a;, we obtain 

/y^ nQn{y) \ (i^Y^ Pn-i{y) \ ^ y^ nPn-i{y) 

n— 1 n— 1 n— 1 

By equating the coefficients of x^"~^, we obtain 

n-l 

(27) nP„_i == nQ„ + ^ kQkPn~k-i, {n > 2). 

fc=i 
Now we substitute the values of the Q„ given in ([25]) 

n-l 

nP„_i = 7iP„ - n{n - l)P„_i + nP,j„i + ^ kQkPn-k-i 

k=l 

so that 

n-l 

nP„ - n^Pn-i - nP;^_, - ^ k{Pk - (fc - l)Pfe-i + P^_i}P„_fc_i. 

fc=i 

D 

From this expression it is very easy to compute the first terms of the expansions 

y-1 y-2 y2_6j/ + ll 2y^ - 21y'^ + 84y - 131 



V = 1 



2x3 6a;4 

3y4 - 462/3 + 294J/2 _ QMy + 1333 



12x5 
y - 1 y^ - 4j/ + 5 2y3 - I5y2 _|_ 42^^ _ 47 



2x2 6x3 

3y4 - 34y3 + i56j/2 _ 366y + 379 



12x4 
Theorem 4.5. M^e /lawe 

(a) For n> I, the degree of Pn is less than or equal to n 

(b) n\Pn{y) has integer coefficients. 

Proof. The equation (j20p may be written 

V-l-^^- (log V-V-xV^~Vy). 

XX 

Since xVic S A, it is clear that 

" Pn{y) 



xV-x-y=-l^Y.-^eA. 



X" 
n—1 



This implies that the degree of P„ is less than or equal to n. 

We prove (b) by induction. The first few P„ satisfy this property. We define 
Pk :— k\Pk so that the recurrence relation (125]) may be written as 



Pn = n^Pn-i - npj^_^ + 



'n — 2 



Hence, by induction, all Pn have integer coefficients. D 
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The most significant contribution by Cipolla is his proof of a recurrence for the 
coefRcients an,k of P„ (see ([30])), which is better than the recurrence given in (f26| . 
We intend to give a shghtly different proof. The result of CipoUa is equivalent to 
the following surprising fact: The solution V of equation (pO)) formally satisfies the 
following linear partial differential equation: 

(28) V^{x-l)Vy-xV^. 

This equation can easily be deduced from the following Theorem. 

Theorem 4.6. For n> 1, we have 

in - i)Pn-i{y) = K-iiy) - Kiy), in > i) 
in~i)Qn-iiy)^Q'n^M-Q'niy). («>2). 

Proof. We will proceed by induction. For n < 3 it can be verified that these 
equalities are satisfied. 

We now assume that (|29p is satisfied for n < N, and we will show that these 
equations are true for n = N + 1. 

By differentiating ([24| with respect to y we get 



/y- Qniy) \ f 1 , y^ Pn-iiy) \ ^ y^ K-iJy) 

Ti— 1 n— 1 n— 1 

so that by equating the coefficients of x^^^^ and of x^^ we obtain 

Af-l N-2 

Qn+1 — Pn ^ /_^ PkQN-k^ Qn — Pn-1 ^ /_^ PkQN~k-l- 

k=0 k=0 

Subtracting these equations we get 

N-2 

Qn+1 ^ Qn = Pn ^ Pn-1 ^ /_^ PkiQN~k ~ QN-k-l) ^ Pn-1 

k=0 

and by the induction hypothesis this is equal to 

N-2 



{N ~ 1)Pn-i +J2Pk-iN-k- l)QN-k-i - Pa 

= -NPn-1 + J2 kQkPN-k^i- 



k=0 

N-1 



fe=l 

By (|27|) this is equal to NP^^i — NQpf so that we obtain 

Qn+1 ~ Qn ~ —NQn. 



This is the second equation of ([29]) for n = iV + 1. In order to achieve the result 
for the first equation, observe that from (P5|) we get 

NQn = NPn - N{N - 1)Pn-i + NPj^_, 
~Q'n = -P'n + iN- 1)P'm-i - P'n-i 
Qn+1 = Pn+1 ^ ^Pn + Pn- 
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By adding these equations we obtain 

- NPn -Pn + P'n+i + ^{^w-i -Pn-{N- 1)^w-i}- 

~ {P'n-i ~P',i-{N- l)P^_,} = NPn -P'm + P'n+i 
which is the first equation of (P^ for n = A^ + 1. D 

We define the coefficients a„ j, imphcitly by 

(30) P„(y) = ^ J, (a„,o;/" - a^ay""' + ' ' ' (-!)"««.„) = 

Analogously, (5„ is of a degree less than or equal to n, and we define the coefficients 
hn,k imphcitly by 

(31) Qn{y)= ^ \ V(-l)H„,fcy"-^ (n>i). 

Remark 4.7. Po{y) has degree 1, which is not given by pop . However, we can 
extend the definition of a{n, k) in such a way that, for n > 1 we have a(n, k) ~ 
for fc < or fc > n. Then a formula such as (|30|) also holds for n == if we add 
up the values from k — —1 to k = n and define a(0, 0) = 1, a(0,— 1) = 1 and 
a(0, fc) = for other values of k 

(30 bis) P„(y)^^^^ ^ (_!)'= a„,y"-'=, (n>0). 

fe=-i 
Note that Qo{y) remains undefined. 

Theorem 4.8. For I < n and < k < n, we have (when defined) 

^(^i ^1) , , ^i(^ ^ 1) , 

(32) an.k = na„^i k-i H ; — an-i,k, bnh — nbn-i,k-i -\ ; — o„_i t. 

n — k n — k 

For 1 < n and < k < n, we have 

(33) On^k = bn,k + {n- k + l)an,k-i- 
For n> 1, we have 

(34) bn.n = nan-l,n-l + /J ( I. ) ^ ^'"''^ "'ri-k-l,n-k-l- 

Proof. ([5^ is obtained by equating the coefficients of y"^-^-^ in the first equation 
in ((29l) . In this way, we obtain 

I'l'^n+k 

(n — 1)! 

^{n- k)— -— — a„_i,fc-i - (n - fc) j a„.fc. 

(n — 1)! n\ 

If n 7^ fc, then the equation for an,k in \Z2\ is obtained. The other equation in bn,k 
is obtained analogously from the second equation in (f29l) . 
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To prove ([33|) . observe that by ([25]), Q„ = Pn — {n ~ l)P„_i + Pn-i, and from 
291) it follows that 



(35) Qn^Pn+P'n, {u > I) . 

Now by equating the coefficient of y"^^^ in both members of this equality we obtain 

Finally (p4)) follows from (l27l) . Recall that — .Sziii and — ^ are respectively the 
values of P„(0), and (5„(0). Hence, by setting y = in ([27)) . we obtain (p4)) through 
multiplication by {n — 1)! and the reordering of the terms. D 



The main problem now is that equations (1321) do not allow us to compute the 
coefficients a„_„. CipoUa gives an algorithm to simultaneously compute the coeffi- 
cients a„.fc and 5„,fc based on Theorem 14.81 In the procedure of Cipolla, these key 
coefficients an^n are recursively computed using all the previous coefficients. We 
prefer a method that computes An := a„_„ and i3„ := ««,«-! separately and then 
compute the remaining coefficients by using p2p . 

Theorem 4.9. In order to compute the numbers a„_fe, we may first compute the 
sequences An '■— On^n o,nd Bn '■— On.n-i by the recursions 

(36) Ao-1, Ai-2, Bo = l, Bi = l, 

(37) Bn = nBn-i+n{n-l)An-i 



n-l ^ „ 

rt — 2 



(38) An = n^An-i + nBn-i- 

n-l 

- (" - 1) E ( r 1 1 ) {^(^ - i)^'«-i - ^'^ + ^^fe-i} A.-fe-i. 

After this one we may obtain a{n,k) := On.k- Setting 

a(0,0) = l, a(0,-l) = l, a(l,0) = l, a(l,l) = 2 
and all other a(0, k) and a(l, k) — 0. Then, for n > 2, put 
a{n,n) = An, 

(39) a{n,k) = na{n-l,k-l) + — -^a(ri-l,fc), (0 < fc < n) 

n — k 

where a(n, fc) = for k < or k > n. 

Finally, we may obtain the b(n, k) :— bn,k from 

(40) 5(71, fc) = a(n, k) - {n - k + l)a{n, k - 1). 

Proof. The constant term of P„ is — ^ and the coefficient of y in P„ is -^, so that 
equation ([57]) follows from the first equation in ([2^ taking it with y = 0. 

In the same way, ([55]) follows from (P5|) . and (15^]) is the first equation in (|32p . 

Equation (HUl) for the b{n,k) follows easily from (|35l) . D 

The array of coefficients a(n, fc) for < n, fc < 7, reads 
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1 























1 


2 




















1 


6 


11 

















2 


21 


84 


131 














6 


92 


588 


1908 


2666 











24 


490 


4380 


22020 


62860 


81534 








20 


3084 


35790 


246480 


1075020 


2823180 


3478014 






720 22428 322224 2838570 16775640 66811920 165838848 196993194 



and the b{n, k) for 1 < n < 7 and < fc < 7 are 



1 


1 




















1 


4 


5 

















2 


15 


42 


47 














6 


68 


312 


732 


758 











24 


370 


2420 


8880 


18820 


18674 








20 


2364 


20370 


103320 


335580 


673140 


654834 






720 17388 187656 1227450 5421360 16485000 32215008 31154346 



Theorem 4.10. (a) The coefficients b(n,k) are integers. 

(b) a{n, fc) > and h{n, k) > 0. 

(c) a{n, fc — 1) < a(n, k) for 1 < k < n. 

(d) Forn> 1, a(n,0) = (n - 1)!. 



Proof, (a) We have proved in Theorem 14.51 that the numbers a{n,k) are integers, 
so that from (HU)) . the coefficients b{n,k) are also integers. 

(b) We proceed by induction on n. Assuming that we have proved that a{m, k) 
and 6(rn, fc) are positive for to < n, it foUows from p2p that a(n, fc) and b{n, k) are 
positive for < fc < n. Then (p4l) imphes that b{n, n) > 0, and (p3| with k — n 
proves that a{n,n) > 0. 

(c) This is a simple consequence of ([55]) . 

(d) The equation follows from (15^ by induction. D 

Theorem 4.11. By means of the rule in Theorem \4.9\ one may compute all coeffi- 
cients Un.k of the polynomials Pn{y) for 1 < n < N in 0{N'^) coefficient operations. 

Proof. We count the operations needed, following the indications in Theorem 14. 9[ 
to compute every a„^fc for < n < iV and < fc < n. 

First we must compute the numbers (™) for < to < iV — 2. Using the scheme 

of the usual triangle, we need to carry out J2k=i ^ additions, which involves {N — 
2)(7V-3)/2 operations. 

The numbers Bn must now be computed for 2 < n < A^ by means of the formula 

Bn^n* {Bn-i + {n - 1) * A„_i). 

Each Bn requires 4 operations, therefore a total of 4(A^— 1) operations are needed. 
We compute the An ioi 2 < n < N using the formula 

An — n * n * An-i + n * Bn-i — [n — 1)* 

* ^ ( , _ j * {fc * (fc - 1) * Ak^i - Ak + k* Bk-i} * An-k-i- 

k=l ^ ^ 
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Hence An requires 7+Y^]^^-^ 8 = 871 — 1 operations. All An together take J2n= 



in- 



1) = AN'^ + 3N — 7 operations. These numbers are the a„^„. The ao^k and ai^k 
require no operations. Finally we compute for < k < n 

an,k = n* {a„_i,fe„i + {n - I) * a„„i^fe/(n - k)}. 
Therefore, a„_fe takes 6 operations. For each n, every Un^k for 1 < k < n takes 
6(n— 1) operations. And each an,k for 2 < n < TV takes X]n=2 6(n— 1) — 3N{N—1). 
The total cost in number of operations is therefore 

^ — ^^ — ^ + 4(iV - 1) + AN'^ + 3N-7 + 3N{N - 1) = 

= ^{15N^+3N-16). 

D 

5. Bounds for the asymptotic expansion. 

5.1. The sequence (a„). First we define a sequence of numbers as the coefficients 
of a formal expansion in A. 

Lemma 5.1. There exists a sequence of integers (a„) such that 

00 , T 00 _, 

(41) log(i-E-) =E--- 

n—l n—1 

The coefficients may be computed by the recursion 

n-l 

(42) ai = 1, On — nl ■ n + y kl On-k- 

fc=i 

Proof. It is clear that u = 1 — J2 nlx^'^ E U C A, so that u^^ E U and logw^^ are 
well defined. To obtain the recursion we differentiate (1411) to obtain 



■^-^ n ■ nl /\-^ a„ \ /_ ^-^ n! 



\l^^n+l)[^ 1^ 



X" 
n—l n—l n—l 



Equation ((42]) is obtained by equating the coefficients of a; " ^. The recurrence 
(|42p proves that a„ is a natural number for each n > 1. D 



The first terms of the sequence (an)5^i are 

1, 5, 25, 137, 841, 5825, 45529, 399713, 3911785, 42302225, . . . 

Lemma 5.2. For each natural number n we have 

(43) a„ < 2n • nl. 

Proof. We may verify this property for ai, 02, 03 and 04 directly. For n > 4 we 
proceed by induction. Assume the inequality for ak with fc < n, so that by (j42p 



rj,' ■ n ^ — ^ 



k /n^ 



nl ■ n ^-^ iji — fc)! • (n — fc) n \k 



< 



k=l 



,fc/ n^ J \n n^ n{n — 1) 

fc=2 ^ ^ ^ ' 
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For n > 4, it is easy to see that this is < 2. 



D 



Lemma 5.3. For each natural number N there is a positive constant cm such that 
(44) a;(i_V:^)>i, x>CN. 



Proof. It is clear that the left-hand side of ((31]) is increasing and tends to +00 when 
X — >■ +00, from which the existence of cjv is clear. 

The value of cjv may be determined as the solution of the equation 

N I 



(45) 



1, 



X > 1. 



In this way we found the following values. 



Cl 


2 


C6 


4.15213 


Cll 


5.61664 


C20 


8.70335 


C2 


2.73205 


cy 


4.43119 


C12 


5.93649 


C30 


12.34925 


C3 


3.20701 


cs 


4.71412 


Cl3 


6.26449 


C40 


16.03475 


C4 


3.56383 


cg 


5.00517 


Ci4 


6.59947 


C50 


19.72833 


C5 


3.86841 


ClO 


5.30597 


Cl5 


6.94035 


C60 


23.42351 



D 



Remark 5.4. Notice that for x > cn the sum in (|44|) is positive and less than 1. 

Proposition 5.5. For each natural number N there exists cIn > such that, for 
a; e C with \x\ > cIn, there exists 6 with \9\ < I such that 

n—1 n—1 



N+lx^+^' 



> dN- 



Proof. By comparing the expansions (j41l) and 

N , , N 



(47) 



it is clear that 



I \ -1 



M-E^)^E 



a„ 1 



E 



&„ 1 



+ (^ + 1)! 



°"^-i^ , so that bj^^i < fljv^i. 



The above expansion is convergent for all sufficiently large \x\, so that 

Ebn ^ _ bN+i 1 , . 

77^ " TV + 1x^+1'^^^^^ 

n=N+l 

where lim^^^oo gN{x) = 1. Hence there exist sufficiently large cIn such that 

\bN+igN{x)\ < ON+l, \x\ > dN- 

This ends the proof of the existence of dN ■ 
We have 



('-(-Ei: 



I \ -1 



n=l 



L. ^ 77 -r 



n=l 



OJV+I 



flAT + l 



E 



„ri-Af-l 
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Since all a„ and 6„ are positive, this is a decreasing function for x — >■ +00, and the 
lowest value of cLn will be the unique solution of 



M'-ES^ 



I \ -1 



N 



1 \{N + l)x^+^ 



n X" 



1. 



a.N+1 



We obtain the following table of values 



di 


1.03922 


de 


4.54145 


dii 


5.73661 


^20 


8.73298 


d2 


2.38568 


dr 


4.75734 


di2 


6.03061 


dso 


12.37349 


ds 


3.33232 


ds 


4.97336 


di3 


6.33969 


^40 


16.05983 


di 


3.92171 


dg 


5.20626 


du 


6.66091 


^50 


19.75448 


d5 


4.28707 


dio 


5.46090 


di5 


6.99175 


deo 


23.45053 



D 



Remark 5.6. The numbers dj^; in Lemma [53] are very similar to the numbers cjv of 
Lemma [5.31 This is no more than an experimental observation, but since the c^ 
numbers are easy to compute and d^ are somewhat elusive, it has been useful to 
start from cn as an approximation to d^ in order to compute d^- 

5.2. Some inequalities. 

Lemma 5.7. For u > 2 we have logali(u) < 21ogM. For u > e^ we have ali(u) < 
2m log w. 

Proof. The first inequality is equivalent to ali(M) < u^. Since li(a;) is strictly in- 
creasing, the inequality is equivalent to u < li(u^). 
For u > 2 we have h(u) > h(2) = 1.04516 ... so that 



\i{u^) = li{u) + 



"" dt 



\ogt 



>1 



logM^ 



Hence, our inequality follows from j" J^ > m — 1, that is from u > 21ogM. However, 
this last inequality is certainly true for u > 2. 

The second inequality is equivalent to u < li(2ulogu) and has a similar easy 
proof. D 



Lemma 5.8. For all integers n > 1 we have 
(48) 



e/- l0g"+4 

where /„ — 4(n + l)/3. 



" (log log i)" (log log u)"- 

dt < iu- 



log"+^ u 



{u> e-^") 



Proof. Notice that /„ > 1. For t > e the function log log t is positive and increasing 
so that 



r (log log i)" , „ , ,„ r dt 

/ \ „+i, dt < log log u)" / -^^ 

Jef^ l0g"+'i Je/„ l0g"+' 



It remains to be shown that 



dt 



< 



4u 



b„n+l I — 1 n+1 ' 
R t log U 



:g"-^ t 



Replacing u by e^ this is equivalent to 



^ e* , 4e^ 
dt< 



i"+i 



„n+l 



(u > e^"). 



(X > In)- 
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For the function 

n('r\ ■= — / 



4e^ 
G(x) ■= —— - / — — dt 






we have 

4(n+ l) 
x"-^^ \ X 

so that for a; > 4(n+ l)/3 we obtain G {x) > 0. Since G'(/„) > we have 0(0;) > 
for aU X > fn- □ 

Theorem 5.9. The polynomials Pn{y) defined in (j2ip satisy the inequalities 

(49) |P„(2;)|<3-n!y", y > 2, n>l 

and \Pa{y)\ <yfory>2. 

Proof. Since Po{y) = y — 1, the second assertion is trivial. 

Given r > 0, for each polynomial P{x) = X!n=o '^n^" ^^ define 



TV 
n=0 



1^11 = Ei«" 



It is easy to show that 

||P + QI1<I1P|1 + I1QI1, I1PQI1<I1P|1-|1QI1 

and that for the derivative of a polynomial of degree < N 

^ AT ^ N 

IIP'II = E n|a„|r"-i < - ^ |a„|r" ^ -||P||. 

n=0 n=0 

For y > r we have the inequality 

Af Af AT 

\P{y)\ = |5] a„y"| < J^ \a,,\y" < J] |a„|r"(y/r)" < {y/rf\\P\\. 

n— n— n— 

Hence, our Theorem follows if it can be shown that for n > 1 we have \\Pn\\ < 
3-2"n! (for r == 2). 

Define 5„ := ||-Pn||- By ()30p we have S*,! = — P„(— 2), and it can be shown that 
S™ <3-2"n! for <n < 15. 

For n > 15 it follows from (|26l) and the aforementioned properties of ||P|| that 

1 """'^ h 

Sn < nSn-1 + -zSn-1 + - ^ fcUfc - l)^*;-! + 5fe + - Sk-l) Sn-k-1 ■ 

k=l 

It follows that Sn < r„ where T„ ~ S„ < 3 • 3" n! for < n < 15, and that for 
n > 15 

3n Iv^V,^. K^k~2) \ 

Tn '■— -^Tn-l H 2^\ ^ ^ 9 Tk-ljTn-k-l- 

fc=l 
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Now we proceed by induction. For n > 15 and assuming that we have proved 
Tfc < 3 • 2*^ fc! for k < n, we obtain 







„ n— 1 

+ -E(fc2^ 


k\ + 


fc(3fc- 
2 


2)2fe- 


-i(fc-l) 


1^2"- 


-»-'(»- 


/c- 


1)! 


Hence 






















Tn 

3-2"n! 


- 4 


3 v-^ / k ■ k\(n 


-k 
■nl 


-'^ 


{3k 


-2)fc!(n 
8-n! 


-k- 


"^). 










3 3 

- 4 ^ 8n2 


>; 
fe=i 


7fc-2 


- 4 


3(7n- 
8n2 


»' + 


o "-2 

3 V . 

8n2 ^ 
fc=i 


7fc- 


-2 



Therefore, by using the symmetry of the combinatorial numbers, we obtain 



n-2 



Tn 3 3(7n-9) 3 ^7n-ll 



<- + i^—^ + -—Y,——i< 



3 • 2" n! - 4 8n2 16^2 ^-^ f""!) 

fc=l \ k J 



3 3(7n - 9) 3 ^, 7n-ll 

- 4 + 8n2 + 16^ ■ ^" " ^^ n-1 " 
3 3(7n-9) 3 ,^ , , n(4n - 63) + 87 , 

for n > 15. D 

Corollary 5.10. We have 

(50) |P„_i(j/)|<n!y", n > 1, y > 2. 

Proof. This foUows easily from the above Theorem. D 

5.3. Main inequalities. To simplify our formulae we introduce some notation. 
First we set r„ := 3 • nl so that, for n > 1, we have |P„(i/)| < r„y" when y > 2. 

Let c„ and d„ be the constants introduced in Lemma 15.31 and Proposition 15.51 
Let a„ be equal to max(e, c„, d„) and let /?« > e be the solution of the equation 

(51) = an- 

log a; 

(The function j-^ is increasing for t > e). 

Finally, define a;„ :— max(/3„, /„, e2), where /„ is defined in Lemma 15.81 

Proposition 5.11. Let x be a real number such that x > Xn, and set y :— logx. 
Then 

y > 2, a; > c„y, x > (i„y, x > fn- 

Proof. Since x > Xn — max(/3„, /„, e'^) we have x > e^ , so that y = logx > 2. 
We also have a; > /?« > e. Since j^^ is an increasing function for t > e we 

obtain j^ > jj|^ = a„ = max(e,c„,d„). Therefore, | > c„ and | > d„ as 
required. D 
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We insert a table of the constants 



Xi 


7.38906 


Xq 


10.81135 


Xll 


16.00000 


a;20 


29.57923 


X2 


7.38906 


X7 


11.70187 


2^12 


17.33333 


2^30 


47.86556 


X3 


7.38906 


Xs 


12.60164 


a;i3 


18.66667 


2^40 


67.69154 


X4 


8.29874 


Xg 


13.58167 


a;i4 


20.00000 


2^50 


88.57644 


X5 


9.77283 


2^10 


14.66667 


2^15 


21.42740 


2^60 


110.29065 



For each natural number N we set 



N 



(52) 



Wn = 1 + J2 



Pn-i{y) 



n=l 



and frequently we write W :— Wn when N is fixed. 

Proposition 5.12. For N > 1 let W = Wn (as in ([5^ ). Then for x > xn and 
y = log a; there exists with \0\ < 1 such that 



(53) 



W + xW + xWx + Wy - X - y - logW = e ■ rN+i 



y 



rN- 



Proof. Denote by T = T{x, y) the value of M^ + xW + xW^ + Wy - x - y -\ogW . 
Then we have 

T = (1 + a;) V '^"~^'^^'' - V '^^""^'^^^ + V I^niM- 

n— n— 1 n— 1 



From ([24)) we have the expansion 

(54) >o8(i+f:^)"^-E%l 



(») 



n=l 



Prom Proposition 15.111 we know that y = log a; > 2 and x > yd^- From (|50|) . for 
y > 2, we have \Pn-i{y)\ < nly" so that we have the majorant 



(55) 



N 



Pn-iiy) 



N 

<log(l-^ 



{ (x/yY 



(by considering this expression as a power series in a; ^, and y as a parameter). 
From dSll and (l55l). we obtain 



N 



(56) 



log 1 + E 



^n(y) 



N 

E 



Qn(2/) 



<S'Ar(x,2/) 



71—1 n—1 

where S]\[{^,y) is a power series majorized by the Taylor expansion of 



Mi-E 

(compare equation (j47| ). 



1 _ Y^ Qn 1 

^1 (a;/2/)"y ^ n (a;/?/)^ 
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By applying Proposition 15.51 we deduce that, for x > ydN, there exists 9 with 
\0\ < 1 and 

(57) ^-(-'^)-^^!^- 

If we substitute ([56|) in the expression for T, then by Theorem 14.31 ^H the terms 
in x^" with n < N cancel out, and the terms in a;~^ add up to — P/v(?/)x^^. It 
follows that 

(58) T = -^^ + SN{x,y). 

x" 

Therefore, since y > 2, we have 






N 



X 

so that from (|43)) , 

x™ V X / x" x" 

where TV > 1 and § + ^^^^ < 3 for a; > e^ are applied. D 



Proposition 5.13. For each natural number N let un — e^" . Then there exists 
vn > UN such that 

(59) HfNiu)) -u = 0- 13{N + 1)! ^(|Qg^Qg")'^ ^ (^ > „^) 

log u 



where \9\ < 1 and 



N 

(60) /Ar(e^) := a;e^W^Ar(a;,logx) = xe^(l + ^ 



F„_i(logj:) 



n=l 



Proof. To simplify the notation, the abbreviation W{x,y) = Wjv(a;,j/) is used. 
Differentiating (|60p we obtain 



^(li(/^(e^))-e-) 



^-——-(e'^W + xe'^W + xe"" (w^ + -Wy) ] - e^. 
Mle^)) I \ X J } 



l0g(/Ar(e-)) 



Assume that x > xat , so that x > d^ log x and x > e^ . We may apply ((53)) to 
obtain 



-(li(Me^)) - e^) = ^^^^^^{ly + .ly + .14^. + M.,} 



This may be simplified to 



'^ hUf (.^W -^-^ - ^"^ ^ T-jv+i log^a; 
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Since x > xn we have x > ycN, so that by Lemma [5?3l 

P„_.i(logx)^ / — "' 



<i + E^^^-f^)|^K^-E 



^1 i'^/yy' 



> 1 



that is xWn{x, \ogx) > 1, so that log(/Ar(e'^)) > x. Hence, for x > xn (with 
another 0), we have 

— (h(/^(e-)) - e) = 0^ e^ 

Defining Hn{u) :— li{fN(u)) — u the above equation is equivalent to 

, rN+i log^a: 
■^^(^ ) = ^ ^JvTi ' (^ - ^^) 

and, since um '■— e^™, 

WW ^ /j'^JV+i (loglogu)^ 

log u 

Lemma 15.81 can be applied since xm > Jn: so that u > un > e-^". Hence, 
integrating over the interval (un,u) we get 

Tj f \ u ( \ . a 4rjv+iM(loglogM)^ 

log u 
The function (A^+1)!t — tJ^^— -(loglogw)^ is increasing (as product of two positive 

increasing functions) for u > e^" , so that there exists v^ > un for which this 
function is greater than Hn{un), so that 

^ 13-(JV + l)!7.(loglogu)^ 
HnW^B , {u>vn). 

log M 

D 

Remark 5.14. For the values of n appearing in our tables, the equality Un — Vn 
holds, since, in these cases, 

„/ (ri + l)!u„(loglogu„)" 

l0g"+' Un 

Lemma 5.15. For any natural number N, and u > e^" we have log fN{u) < 
2\ogu. 

Proof. First observe that the hypothesis u > e^" implies (with u — e^) that x > 
Xn, so that logx > 2 and x > cjm logx. fProposition lS.lip . 

The inequality log /at (u) < 21ogu is equivalent to fN{u) < u'^, and together 
with u = e^ it is equivalent to 

Pn~lil0gx)\ ^ 2:e 



-(l+E 



xe-^ I 1 + > ^^^^^^ I < e^ 



From Corollarv l5.10[ since x > e^ and x > cat?/, and by Remark l5.14[ 
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Hence our inequality follows from 

X N I 

ye X ( v-^ n\ \ 

Since we assume that x > ycN, the second factor is greater than 1, so that 

X y V _ {x/y)^ J X 

Finally, it is easy to prove that e^ log a; > x^ for x > e^. D 

The asymptotic expansion with bounds can now be proved. 
Theorem 5.16. For each integer N > 1 

(61) ^l[(u)^fj,{u) + 260iN + iy.u(^-^f^^^^Y, {u>vm), 

\ log u / 

where vn is the number defined in Proposition \5.13\ . 

Proof. Since li(ali(u)) = u, Proposition 15 . 1 3l vields. for u > vn, 

liiMu)) - h(an(.)) = / -£- = miN + 1)! u^ip^l?^. 

Jali(«) fog I log ^ U 

Since vn > un = e^", u > vn implies logw > 2, hence u>2. 

From Lemma 15.71 logali(w) < 21ogw, for u > 2. Analogously, Lemma 15.151 
implies that log /at (m) < 21ogu, for u > e^™. Therefore, for u > ujv, we have 



|ali(u)-/Ar(u)| ^ 



/iv(u) ^^ 
ali(«) log* 



21ogu 

It follows that there exists 9' with \9'\ < 1 such that 

r/iv(u) ^^ 



ali(«)-/jv(u) = 0'(21ogu) 

Ja,li(u) ^Ogt 

and the result follows easily. D 

The actual error appears to be much smaller than that given in Theorem 15.161 
However, as usual with asymptotic expansions, having a true bound allows realistic 
bounds to be given of the remainder for specific values of N. 

The true error after N terms of an asymptotic expansion, while the terms are 
decreasing in magnitude, is often of the size of the first omitted term. In our case, 
the magnitude of the term P/v(logx)a;~^~^ depends on the polynomial PN{logx). 

Numerically, it appears that for n > 3: 

(62) \Pn{y)\<(-, ) y", (y>21ogn) 

V e log n / 

although we have not been able to prove this. 

From Theorem 15.161 more realistic bounds can be obtained for the first values 
of N. This is done in the following Theorem. 
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Theorem 5.17. For 2 < N <11, we have 

n=l 



(^^) ^^ = 1 + 2^ -n +^-20-i^to^J •^^^' ^">"^) 

where 



Z2 = 1.50, Z3 = 2.34, Z4 = 3.32, zg = 4.33, zg = 5.36, 
Z7 = 6.39, zg = 7.43, zg = 8.46, zio = 9.50, zn = 10.53. 

Proof. By taking A^ = 10 in Theorein l5.16[ we have, for u — e^ > e^'-", (recall also 
Remark ??) 

ali(e^) -^ P„_i(loga:) log^°x 

^^ = 1 + 2^ ^ + ^^^T^' (^>2:io) 

n—1 

with i? = 26 • 11! = 1 037 836 800. 

We compute the maximurr|3 A/„ of |P„_i(loga;)/log"^^ a;| for a; > xio, so that 
for any 2 < A^ < 10, we have 

ali(e"^) ^ -^ ^ y^ P„-i(loga:) ^ 



log^/ f. P„_i(log:r)log"-^^-^x 1 



n=Af+l 



SO that 




aH(e^) 


^ 1 y> P„-i(loga;) ^ 




, ,log~x^ ^ Af„log"-^-^x ^ „W"-^x 




n=Ar+i 



We determine a value z^ > xio such that, for x ~ z^, 

n=Ar+l ° 

Since this is a decreasing function of a;, we obtain for x > z^ 

(64) «:l = l + V^i.iMfi)+..20(-^)"i^. 

n— 1 ^ 

We consider the function 

N 



ali(e-) ^_y.P„-i(logx)^ a;^+i 



on the interval (1.3, z^), to determine the least value of zm for which (|M|) is true. 



3a/2 = 1, A//3 = 1/2, M4 = 1/3, Afs = 0.250636, Afg = 0.526887, Afy = 1.300565, Afg = 
3.719653, Afg = 12.070813, Afio = 43.788782. This last maximum would be much smaller if the 
maximum were taken from a point slighthly greater than xiq. 
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In this way we find: Z2 = 32 and then Z2 = 1.5; Z3 = 49.5 and then Z3 = 2.3395; 
Z4 = 82 and then Z4 — 3.3114; Zg — 155 and then Z5 = 4.3237. 

If we take TV = 20 in Theorem WM. we obtain z^ = 113, z'-j = 143, Zg = 187, 
z^ = 251 z[q = 353, z^i = 528 from which zg = 5.3514, Z7 = 6.3851, zs = 7.4208, 
zg = 8.4566, zio = 9.4914 and zn = 10.5251 are obtained. D 

Remark 5.18. We have proved ((63|) only for 2 < iV < 11, ahhough something 
similar appears to be true for the general case. If (|63p were true for all n, then for a 
large u = e^ we could take N fn x terms in the expansion and in this way the error 
would be « -2%, so that ali(u) could be computed with an error less than sa 20. 

In fact, for several values of u, the terms of the expansion have been computed 
up to the point where these terms start to increase. Always the computation is 
terminated when N ^ x and the error appears to be bounded. ( For example, 
with u — 10^'^*', we compute 230 terms of the expansion, which coincides with 
loglO^'^'^ « 230.259. The approximate value obtained for ali(u) has an absolute 
error equal to 40.94738, which can be compared with the fact that ali(10^°°) has 
103 digits ). 

6. Applications to p„. 

6.1. Asymptotic expansion of p„. Inequalities for the n-th prime number can 
be found in [18], [16], [12], [4]. In fact, from 7r(a:) = li(a;) + 0{r{x)) we may obtain 
Pn = ali(n) + 0(r(n log n) log n), if r{x)/x is sufficiently small. For example, in 
[l2] , it is noticed that from a result of Massias [10] , it follows that 

(65) Pn = ah(n) + Oine-"^^"^) 

so that the asymptotic expansion of ali(n) is also an asymptotic expansion for p„, 
that is, 



N 



Pfc_i(loglogn)\ / /log log nN^ 



V ^ log^n ) V V logn 

By assuming the Riemann hypothesis, Schoenfeld [H] has proved 

(67) |7r(a;) - X\{x)\ < — V^logx, {x > 2657). 

87r 

This result will be used to obtain (under RH) some precise bounds for p„ 

Lemma 6.1. We have -^^(loga::)^ < ali(x) for x > 94. 

Proof. The inequality is equivalent to 

h(^/a;(loga;)5) < x. 

Differentiating the function f{x) :— x — li(-/i(loga;)5) we get 

^v^^_l 1 ^ log'^'^ , 51og^/^ x^ 

log(Vi(logx)i)V 2^ 2Vi 

Hence this derivative is positive if and only if 

, 3/2 /logx 5\ ^/loga; 5, , 
log3/2a;|^_|_ + _j <7^(^_|_ + _ioglogx 

For a; > 94 we have (log x) 2 < y/x, so that /' {x) > for x > 94. 

Finally, one may verify that /(94) > 0. Hence f{x)>0 for a; > 94. D 
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Theorem 6.2. The Riemann hypothesis is equivalent to the assertion 

1 5 

(68) \pn — ali(n)| < —\/nlog'^ n for all n > 11. 

TT 

Proof. First we assume the Riemann Hypothesis and prove (|68|) . Let r(x) := 
g^^Y^loga;, f{x) := \i{x) — r{x), and g{x) := h(a;) + r{x). For x > 1 we have 
f{x) < \i{x) < g{x), where the three functions are strictly increasing. From (j67p 
for X > 2657, we also have f{x) < tt{x) < g{x). 

The inverse functions satisfy g^^{y) < ali(y) < f^^{y), and if y = n > 7r(2657) — 
384 is a natural number, then g~^{n) < p„ < f^^{n). It follows that the distance 
from ali(n) to p„ is bounded by 

\pn — ali(n)| < max(/~"^(n) — ali(n), ali(n) — g~^{n)). 

Hence, we have to bound f~^{y) — ah(2/) and ali(y) — g~^{y). 
We consider y as a parameter and set a — ali(y), so that li(Q;) — \\{a\\{y)) — y. 
Consider the function m(^) := /(^) — li(a) = f{0^ Vi which is strictly increasing 
and satisfies 

uiO=im-r{0-Ho.)= (^ ^-r{i). 



Therefore, u[a) — —r{a) < and 

u{f-\y)) - f{f-\y)) - li(a) = y - h(ah(y)) = 0. 

If apoint b is found where u{b) > 0, then a < f^^{y) < b, so that b—a > f^^{y)~a 
and one of the required bounds is obtained. 

Therefore, we try b = a + c^{\ogy)2 with c < 1. We have 

u{a + Cy/y {log y)^ ) ^ / j — - - r{a + Cy/y{\ogy)^) > 

> , , ^7^ ;t- - r{a + cV^(log y)^). 

log{a + Cy/y{\og y) 2 ) 

From Lemma |6T] for y > 94, we have y/y{logy)^ < ali(j/) — a, so that 



^a + cVy(logy)i)>^^f|^-K2a) 



c^i\ogy)-2 ~ ^V2^\og\2a) 



log(2a) 

We want to show that this expression is positive. For y > 94, we have a — ali(y) < 
2y\ogy (by Lemma [??7|) . so that a < 2y\ogy < Ay logy < y^ (for y > 94), which 
yields (with c = l/vr) 

5 I , — 2 
Cy/y{logy)^ - — V2alog (2a) > 

OTT 

> Cy/y {log y)^ - — y^4y log y log^ (4?/ log y) > 
ovr 

> cVy(log2/)^ - — V4ylogylog^(2/^)) = 0. 

ovr 

Hence, we have proved that f^^{y) — a < ^ ^yy (log yY^^'^ for y > 94. 
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To bound a — g^^{y), we consider the function v{^) :— g(^) — li(a) = 5(0 ^ V- 
Then 

and v{a) ~ r{a) > 0, v{g^^{y)) — g{g^^{y)) — y — 0. If a value b is found such 
that v{b) < 0, it wih follow that a — g~^{y) < a — b. 

Choose 6 = a — Cy^(logy)^/^ with c = i. We claim that v{b) < 0. We have 

v{b)=r{b)- r-^<r{b)~^ 
Jb logi log a 

and our claim will follow from r(b) log a — c-^y (log y)^^^ < 0. Finally, since b < a = 
ali(y) < 2ylogy, and by Lemma ISTZl 

1 r^T^ o 1 



r(5)loga < —y/ali{y)(logal[{y)f < —y/2y\ogy{2\ogyf 

OTT OTT 

which proves our claim. 

Hence, (by assuming RH), we have proved that |p„ — ali(n)| < i-y/n(logn)2 for 
n > 384 > 94. By verifying all 1 < ti < 385, we find that the inequality holds 
except for n < 11. 

The reverse implication is simple. From |p„ — ali(ri)| = 0{n^^^) for any e > 0, 
we may derive that 7r(x) = li(x) + 0{x^^^). It is well known that this is equivalent 
to the Riemann Hypothesis. D 



Remark 6.3. The inequality (j68p is only proved by assuming the Riemann Hypoth- 
esis, but is stronger than those contained in [18], [16], [12], [4]. Inequality (|68l) gives 
approximately half of the digits of p„. If our conjecture that ([63]) is true for all N 
is also assumed, then the asymptotic expansion gives about half of the digits of p„. 

6.2. Inequalities for the n-th prime. Let 

f \ 1 /i , Y^ Pfe-i(loglogn)\ 
sn[ti) = nlogn 1 + > — v 

where sq — nlogn. 

Cipolla noted that for fc > 1, Pu{y) = (-l)''+^%- + - • • , andPo(2;) =J/-1- Hence, 
except for the first term, eventually the sign of the fc-th term Pfe_i (log log n) log^ n 
becomes (—1)'^. The asymptotic expansion implies that there exist rj.^ such that 

p„ > so(n), n > ro, p„ > si(n), n > ri, 

Pn<S2N[n), n>r2N, Pn> S2N+l[n), n>r2N+l- 

In fact, ro = 2 is the main result in [IS], ri = 2 is proved in [4] and r2 = 688 383 is 
proved in [5]. The value of r^ ioi N > 3 has not been determined. See Theorem 
16.41 for an estimation of r^ by assuming RH. 

The above reasoning may give the impression that the terms of the asymptotic 
expansion of ali(w,) are alternating in sign, starting from the second term. However 
this is not true. For example, computing the first 230 terms for ali(10^°°), we found 
only three positive terms Pq/x, Pi/x^, and P^/x'^. In fact, the sign of the fc-th term 
is that of Pfc-i(loglogn). Thus we are interested in the sign of these polynomials. 
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The polynomials P/v(j/), for 1 < iV < 23 of odd index, have one and only one 
real root, which is positive. Starting from Pi{y) which vanishes at j/ = 2, these 
roots are 

2, 4.23415, 5.83131, 7.43591, 9.07979, 10.6881, 12.2538, 
13.7876, 15.2977, 16.79, 18.2683, 19.7353. 

The polynomials P2, P4 and Pq have no real roots, and all Pg, . . . , P22 have two 
positive real roots. These pairs of roots are: 

(6.4306,8.2185), (7.16158,9.88528), (7.90293,11.4752), 

(8.63359,13.0241), (9.3507,14.5452), (10.055,16.0458), 

(10.7478, 17.5307), (11.4307, 19.003). 

For example, P9(loglogn) is positive only for n > exp(e^''^'), which is a very 
big number. 

The even terms at first sight appear negative. However Pio(loglogn), for exam- 
ple, is negative except in the interval exp(e^'^^') < n < exp(e^'^*'). 

In a certain sense, the inequalities (|69| are the wrong inequalities. These in- 
equalities would hold only for very large values of r^, especially when we want a 
lower bound of p„ > S2N+1 (except for the three known cases). We estimate r^. 

Theorem 6.4. Let r^ be the smallest number such that 

log log n — 2 



(70) p„ > S3 := 7ilogn + ?i(loglogri - 1 



logn 

(loglogn)^ — 61oglogn + 11 
21og^n ' 



n > r^. 



Then, if the Riemann Hypothesis is assumed, 

39 X 10^9 < ra < 39.58 x 10^^ 
Proof. By Theorem 16.21 there exists 9i with 16*11 < 1 such that 

p„ = ali(n) -1- 6*1 -^- — (log 71)2, n > 11. 

TT 

By Theorem I5.17[ with 5 < iV < 10 and setting n — e^, x — logn, and y 
log log n, we have for n > e^™ 

y-1 y-2 y2_6y + ll 



aU(n)=a;e^fl + ^^ + 

V X 



X X 2x^ 



N 



Psiy) , ^Pfc_i(y) y 

1 h f 2CJV 



E 



N 



fc=5 

The inequality of the Theorem is obtained if 

V X^ ^-^ X'^ X"^'-/ TT 

k=5 

This is equivalent to 

P3(2/) 4-5:Pfc-i(2/)e-('=-4)^ > c^y^e-(^-3)^ + ie^^e"^^". 

k=5 
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Since ^3(2/) — > +00 and all the other terms tend to as y — > +cxi, it is clear that 
the inequality is true for y > yo. With N = 10, we find yo = 4.254946453 . . . The 
inequality p„ < s„ is true for n > 3.95702224148845656 x 10^°. This proves that 
rg < 39.58 X lO^^. 

In order to show that rg > 39 x 10^^, we directly show that, for n = 39 x 10^^, 
the opposite inequality p„ < S3 is obtained. 

We compute 

S3 = 2.87527 18639 02974 79681 42399 35057 89294 02005 87915 X lO^^ 

Now we can compute ali(n), for which we already have obtained a good approxi- 
mation through the asymptotic expansion, and then apply the Newton method 

ali(n) =2.87527 18639 02495 21516 14800 14732 45414 39731 X 10^^ 

Therefore, from Theorem 16. 2[ we obtain p„ < ali(n) + i^/nlog^ n, so that 

Pn < 2.87527 18639 02756 97808 39055 05640 30082 86370 11482 X 10^^ 

and we can conclude that p„ < S3. D 
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